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Descripti n 

TECHNICAL FIELD OF THE INVENTION 

s The present invention is directed, in general, to mass flow rate meters incorporating measurement of Coriolis 

forces and, more specifically, to such meters that, by virtue of their geometry, associated compensation circuitry or 
both, are rendered substantially insensitive to fluid pressure and density 

BACKGROUND OF THE INVENTION 

10 

During the last 15 years a new flow meter technology has evolved using vibrating conduits and Coriolis forces 
developed in fluids flowing therein to measure the mass flow rate of a fluid. A typical sensor employs two parallel 
conduits usually bent into some unique shape such as a "U" or an "S" or a bowtie shape. These conduits are normally 
welded to a rigid casting at their inlets and outlets and the conduits are forced to vibrate with one conduit vibrating in 

is opposition to the other for balance purposes. Fluid flowing through the vibrating conduits therefore experiences this 
vibrating motion and, in particular, the angular rotation of the conduits, causing Coriolis forces to impinge on the walls 
of these conduits. The net effect is a slight deformation and deflection of the conduit proportional to the mass flow rate 
of the fluid, the angular velocity of the conduits and some conduit geometry constants. This deformation exemplifies 
itself and is normally measured as a small phase or time delay between the deflection at the inlet ends of the conduits 

20 compared to the deflection at the outlet ends. 

Since its inception and reduction to practice in industry, many innovative methods and inventions have been em- 
ployed to maximize these tiny deflections and their time delays to improve the sensitivity of these devices. Many of 
these efforts have focused on improving the shape of the bent conduits to maximize angular rotational velocity while 
minimizing the stiffness that opposes the Coriolis deflections. This line of reasoning has led to many innovative shapes 

25 that lengthen these phase delays, however the design tradeoffs are generally, (a) a complex conduit shape to bend, 
(b) thinner walled flow conduits, (c) a lower frequency of operation, creating susceptibility to ambient vibrations and 
(d) a higher pressure drop through the device due to introduction of flow splitters and complex conduit geometry. Some 
examples of these flow conduit geometries are the dual S-conduits of U.S. Patent numbers 4,798,091 and 4,776,220 
to Lew, the ^-shaped conduits of U.S. Patent number 4,852,41 0 to Corwon, et al., the B-shapes conduits of U.S. Patent 

30 number 4,891 ,991 to Mattar, et al., the helically-wound flow conduits of U.S. Patent number 4,756,198 to Levien, figure- 
8 shaped flow conduits of U.S. Patent number 4,71 6,771 to Kane and the dual straight conduits of U.S. Patent number 
4,680,974 to Simonsen, et al.. These complex conduit shapes also contribute to the overall size of the device and 
usually preclude the ability to place the device within a reasonably-sized pressure-containing case for safety reasons. 
In addition to the problems associated with the general shape of the conduits, the magnitude and placement of masses 

35 on the conduits (such as magnets), for the purposes of driving and sensing the requisite vibratory motions, can ad- 
versely affect the response of the device to flow rate under the conditions of changing pressure, density or both. 

Currently, several new types of Coriolis flow meters are being developed using radial vibratory motion of the wall 
of a single straight conduit, with flow going either within or around the conduit, thereby eliminating the flow splitters 
and the second counterbalancing flow conduit. Such meters are referred to as "radial mode meters." 

40 a Coriolis mass flow meter sensor is a device that accommodates the flow of fluid within or around one or more 

flow conduits and subjects that flowing fluid to an angular rotation by virtue of a vibration or oscillatory motion of the 
flow conduit(s). The flow conduit(s) are normally designed to be highly resonant devices that can be easily excited to 
vibrate in one or more natural modes of vibration for this purpose. This angular rotation imparted into the flowing fluid 
by this conduit motion thus causes resultant Coriolis forces to develop and bear against the sidewalls of the flow conduit 

45 in a direction that is 90° from both the angular rotation vector and the fluid velocity vector. These resulting forces are 
cyclic in nature (normally sinusoidal) at the frequency of the excitation motion. 

From WO-A-9214123 and US-A-5 373 745 a flow meter apparatus for measuring the mass flow rate of a fluid 
using the Coriolis principal is known. The known apparatus uses a single straight flow conduit which is vibrated in a 
radial mode. Sensors are provided on the flow conduit in order to measure alternate vibrations. The drivers and sensors 

50 may be advantageously located along the conduit such that sensitivity to mechanical pressure or temperature variations 
is minimized. 

The US-A-5054326 discloses a further known Coriolis meter, wherein the conduits having a curved form is vibrated 
in two modes of oscillation for compensating density. 

GB-A 2 171 200 discloses a mass flow meter on Coriolis principal using two parallel straight measuring tubes 
55 wherein the mechanical stresses can be corrected. 

The US-A-47561 98 discloses a further mass flow meter making use of Coriolis effects, wherein a specific mechan- 
ical structure is provided for having specific density sensitivity. 

US-A-5275061 describes a Coriolis mass flow meter in which velocity or acceleration sensors are mounted in 
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special positions to avoid errors due to mechanical, pressure or temperature influences. 

Since these devices are normally designed to be highly resonant and compliant structures, the resulting sinusoidal 
Coriolis force distribution results in a sinusoidal (as a function of time) deflection of the structure that is proportional to 
the mass flow rate of the fluid and the frequency response of the structure for that given excitation. As stated previously, 

s the frequency response of the device can change as a function of fluid pressure, density, temperature and conduit 
stress, it is essential that it be accurately known and compensated for, or held to be constant to achieve accurate mass 
flow rate measurements from the device. Therefore, one aspect of the present invention is the ability to control the 
response of a flow conduit geometry to achieve insensitivity to these changing parameters by virtue of the flow conduit's 
geometry in combination with a motion measurement method. This arrangement implemented on a radial mode meter 

io has many advantages over the bending mode meter types previously described and virtually solves the design tradeoffs 
(a) through (d) previously mentioned. While both radial mode and bending mode meters can exhibit pressure sensitivity, 
the problem is generally more severe on radial mode Coriolis meters because internal fluid pressure tends to stiffen 
the walls of the conduit with regard to radial vibrations more than with regard to bending mode vibrations. In addition, 
radial mode meters generally have vibration frequencies that are orders of magnitude higher and have flow-related 

'5 phase or time delays that are many times smaller than traditional bending mode Coriolis flow meters. For example, a 
typical 2.5 cm (1 ') bending mode Coriolis flow meter operates at approximately 100 Hz while a typical 1 " radial mode 
meter operates above 3000 Hz. Consequently, the Coriolis signal-processing method of choice traditionally has been 
a phase or time delay measurement. Unfortunately, measuring phase or time delay results in very poor sensitivity and 
resolution for radial mode meters. For example, a typical time delay measured at a nominal flow rate on a traditional 

20 bending mode flow meter is in the tens of microseconds and easily measured with today's electronics. By contrast, the 
same flow rate in a comparable radial mode Coriolis meter might cause only tens of nanoseconds of time delay, a 
period of time that is much more difficult and expensive to measure accurately. To date, the only radial mode Coriolis 
flow meter patent has been that of Lang (U.S. Patent number 4,949,583) who describes a radial mode vibrating sensor 
However, Lang exclusively uses a relative phase displacement measuring technique requiring signal processing cir- 

25 cuitry designed to discern minuscule time delays. Such circuitry is relatively expensive and inaccurate, damaging the 
efficacy of the meter as a whole. Consequently, to date no practical device has yet been disclosed or marketed for 
commercial applications. 

Therefore, the current invention addresses the more general task of providing a Coriolis flow meter that is insen- 
sitive to both density and pressure at the same time. 

30 The present invention, which is defined in claims 1 and 15 hereinafter, was therefore discovered, developed and 

applied to solve the pressure, density and small signal problems with regard to prior art radial mode Coriolis flow meters 
and was subsequently found to be relevant to traditional bending mode meters as well. 

Using the present invention therefore, a Coriolis mass flow meter is provided that solves the previously-mentioned 
design tradeoffs and is insensitive to changing fluid parameters of pressure, density, temperature, or any combination 

os thereof by using either of three embodiments to achieve this end. 

Briefly, the present invention aims to provide a Coriolis mass flow meter that is made to be insensitive to pressure, 
density, conduit stress or temperature, by virtue of its geometry and motion detection method or its geometry and 
motion detection method in combination with a mathematical algorithm or its geometry and motion detection method 
in combination with direct calculation of and compensation for the pressure and density of the fluid. 

*o Several embodiments are disclosed that incorporate the desired insensitivity properties but vary in the complexity 

of either the sensor (forming a mechanical portion of the device) or the electronics (forming a signal processing portion 
of the device), since there is a natural tradeoff between the two. 

As previously mentioned, the resulting flow-related deflection is usually a sinusoidal function of time, thereby caus- 
ing a velocity, acceleration and a time delay also to be associated with the motion, all of which can be measured using 

45 fundamental sensing techniques (accelerometers, velocity sensors, displacement sensors, etc.) . As a fortuitous con- 
sequence of nature, the measured response of a given conduit geometry to changes in fluid parameters (pressure, 
density, etc.) varies depending on which motion sensing technique is used. Furthermore, since there is a definite 
mathematical relationship between acceleration, velocity, displacement and time delay, a continuous characteristic 
curve can be defined describing the change in the response of the device to pressure, density, etc., as a function of 

50 sensing technique. The shape and zero intercepts of these characteristic curves for changes in pressure, density, 
temperature and conduit stress, can be effectively controlled according to the present invention, allowing the designer 
to choose the desired response to the extent that a sensor can be designed that is inherently insensitive to changes 
in these parameters. 

In addition, an algorithm can be derived from these characteristic curves that can be advantageously applied to 
55 allow the designer to work at any desired point (defined for purposes of the present invention as a "working point") 
along the response curve. This aspect of the invention allows the designer to use sensor geometries that are not 
inherently insensitive to changes in fluid parameters and by applying the appropriate algorithm of the present invention, 
render its resulting flow-related signals insensitive to changes in fluid parameters. 
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Another embodiment of the invention is disclosed whereby a second mode of vibration is excited to determine its 
natural frequency. Due to a fortuitous relationship between pressure and density and the frequency of these two modes 
of vibration, a simultaneous solution of both the pressure and the density of the fluid can be attained. Knowing the 
value of both pressure and density in conjunction with the characteristic curves just described, a compensation value 
can be directly determined and accurately applied thereby rendering the resulting flow-related signals to be insensitive 
to changes in fluid parameters. 

The foregoing has outlined rather broadly the features and technical advantages of the present invention so that 
those skilled in the art may better understand the detailed description of the invention that follows. Additional features 
and advantages of the invention will be described hereinafter that form the subject of the claims of the invention. Those 
skilled in the art should appreciate that they may readily use the conception and the specific embodiment disclosed 
as a basis for modifying or designing other structures for carrying out the same purposes of the present invention. 
Those skilled in the art should also realize that such equivalent constructions do not depart from the scope of the 
invention as set forth in the appended claims. 

BRIEF DESCRIPTION OF THE DRAWINGS 

For a more complete understanding of the present invention and the advantages thereof, reference is now made 
to the following descriptions taken in conjunction with the accompanying drawings, in which: 

FIGURE 1 illustrates a cutaway side view of a mechanical portion of a Coriolis mass flow meter device of the 
present invention; 

FIGURE 2 illustrates a side view of a vibrating portion of the Coriolis mass flow meter device of FIGURE 1 shown 
in its undetected position; 

FIGURE 3 illustrates a side view of a vibrating portion of the Coriolis mass flow meter device of FIGURE 1 shown 
in a deflected position representing a first elliptical radial mode; 

FIGURE 4 illustrates a side view of the vibrating portion of the Coriolis mass flow meter device of FIGURE 1 shown 
in a deflected position representing a second elliptical radial mode; 

FIGURE 5 illustrates a side view of the vibrating portion of the Coriolis mass flow meter device of Fl GURE 1 shown 
in a deflected position representing a third elliptical radial mode; 

FIGURE 6 illustrates a side view of the vibrating portion of the Coriolis mass flow meter device of FIGURE 1 shown 
in a deflected position representing a fourth elliptical radial mode; 

FIGURE 7 illustrates a graph of curves representing a change in sensitivity of the Coriolis mass flow meter device 
of FIGURE 1 to changes in pressure, temperature, axial stress and density, where some of the curves intersect 
each other at a zero value and at a working point corresponding to a fundamental measurement method; 
FIGURE B illustrates a graph of curves representing a change in sensitivity of the Coriolis mass flow meter device 
of FIGURE 1 to changes in pressure, temperature, axial stress and density, where some of the curves intersect 
each other at a zero value, but not at a working point corresponding to a fundamental measurement method; 
FIGURE 9 illustrates an exemplary circuit diagram for processing the signals of the present invention; 
FIGURE 10 illustrates a graph of actual data taken from a prototype Coriolis mass flow meter device of the present 
invention depicting a percentage deviation in accuracy of the device due to changing pressure; 
FIGURE 11 illustrates a graph of actual data taken from a prototype Coriolis mass flow meter device of the present 
invention depicting the accuracy of the resulting mass flow rate calculations after using the algorithm described 
herein to compensate for pressure and density; 

FIGURE 1 2 illustrates a graph of curves representing a change in sensitivity of the Coriolis mass flow meter device 
of FIGURE 1 to changes in pressure and density, where both curves intersect each other at a non-zero value; 
FIGURE 1 3 illustrates a graph of curves representing a change in sensitivity of the Coriolis mass flow meter device 
of FIGURE 1 due to changes in pressure and density, where the curves intersect each other at a non-zero value; 
FIGURE 14 illustrates an isometric view of two mathematical surfaces describing a range of natural frequencies 
for the first and third elliptical radial modes of vibration, resulting from various combinations of pressure and density; 
FIGURE 15 illustrates a view down a frequency axis of FIGURE 14 showing a relationship and an intersection of 
the two constant frequency lines for the first and third elliptical radial modes of vibration as functions of pressure 
and density; 

FIGURE 16 illustrates an exemplary circuit diagram for processing the signals of the present invention; 
FIGURE 17 illustrates a portion of an exemplary circuit diagram for processing the signals of the present invention; 
FIGURE 18 depicts pressure and density response curves for several Coriolis mass flow meter conduit lengths 
without added masses; 

FIGURE 1 9 depicts pressure and density response curves for the conduits of FIGURE 18 with added center mass- 
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FIGU RE 20 depicts pressure and density response curves for the conduits of FIGURE 1 9 with added outer masses; 
FIGURE 21 depicts pressure and density response curves for a conduit operating in a first bending mode of vi- 
bration without added masses; and 

FIGURE 22 depicts pressure and density response curves for a conduit operating in the first bending mode of 
5 vibration with added masses. 

DETAILED DESCRIPTION 

FIGURE 1 is illustrative of the mechanical portion (a sensor) of the present invention and will be used to describe 
io several embodiment of the invention. The sensor consists of a single straight flow conduit 1 that is made of a strong 
resilient material such as 316 stainless steel or titanium. The flow conduit 1 is fixedly attached at both ends to manifolds 
4, 5 by a welding or brazing process so that fluid flow 21 can easily pass into, through and out of the flow conduit 1 
without leaking. Mounted In association with the flow conduit 1 is a support bracket 2 that is used to mount coils of 
magnet/coil pairs 10 through 17, support any requisite wiring and convey pipeline stress around the flow conduit 1, 
is among other purposes. The support bracket 2 is preferably made of a strong material such as steel pipe and is preferably 
attached at its ends to the manifolds 4, 5 by a welding or brazing process. 

Mounted in association with the support bracket 2 and the flow conduit 1 are magnet/coil pairs 10 through 17. The 
magnets of the magnet/coil pairs 10 through 17 are fixedly attached to the flow conduit 1° by a brazing process, although 
any suitable attachment method can be employed, such as adhesives or fasteners. The coils of the magnet/coil pairs 
20 10 through 17 are fixedly attached to the support bracket 2 by suitable fastening hardware. The magnet/coil pairs 10 
through 15 are used as motion detectors and being magnets and coils, will detect the velocity of relative motion between 
each magnet and its associated coil. This velocity will appear as a voltage on the coil wires proportional in magnitude 
to this relative velocity. The magnet/coil pairs 16, 17 are used as drivers to excite the flow conduit 1 to vibrate in a 
prescribed mode of vibration. 

25 Also mounted in association with the flow conduit 1 are a strain gage 1 9 and a temperature sensor 1 8. The strain 

gage 19 is preferably a semiconductor type or metal film type strain gage mounted on a polyimide backing material 
and arranged to sense the axial (in the direction of flow 21 ) strain in the flow conduit 1 due to temperature or pipeline 
induced axial strains. The temperature sensor 18 is preferably a 100n platinum resistance thermal device ("RTD"), 
mounted to the flow conduit 1 by means of thermally conductive adhesive and arranged to accurately sense the tem- 

30 perature of the flow conduit 1 . 

Concentrically mounted around the support bracket 2 and the flow conduit 1 is a pressure tight case 3 that is 
preferably fixedly attached at its ends by welding or brazing to the manifolds 4, 5 so that an annular space 8 is enclosed 
between the flow conduit 1 and the case 3, that can be held at a vacuum or a predetermined amount of pressure. The 
case 3 is preferably round and made of a strong material such as steel pipe or tubing capable of withstanding internal 

35 pressures or vacuums and capable of withstanding pipeline forces and stresses. 

Mounted in association with the case 3 is a feed through 20 that is used to convey signals through the pressure 
tight case 3 while maintaining a vacuum or predetermined pressure in the annular space 8. Preferably, the feed through 
20 is made of steel, is attached to the case 3 by a welding or brazing process and incorporates glass to metal seals 
for each individual signal carrier to maintain the annular space 8 at a vacuum or predetermined pressure without leaking. 

40 FIGURE 1 , as just described, is therefore a representative Coriolis mass flow sensor geometry that will be used 

to further describe the aforementioned characteristic response curves to changes in pressure, density, temperature 
and conduit stress, the control of these characteristic response curves to achieve insensitivity to these parameters and 
the mathematical algorithm derived from these curves that can be advantageously applied to effect insensitivity to 
these parameters. It is important to note, however, that the methods described herein to control these characteristic 

45 response curves and the algorithms later defined, apply to a wide variety of Coriolis mass flow sensor geometries 
including traditional bending mode meters, as well. The geometry of FIGURE 1 is, however, the preferred sensor 
geometry. 

FIGURE 9 is representative of an electronic circuit that can be used in conjunction with the sensor of FIGURE 1 
to excite the requisite oscillatory vibration on the flow conduit 1 and to extract from the resulting motion, signals that 

so are proportional to the mass flow rate of the fluid flow 21 . 

Since the sensor of FIGURE 1 is designed to be a highly resonant structure, it will normally have the ability to 
vibrate in several natural modes of vibration. FIGURES 2 through 6 depict a progression of natural mode shapes that 
can be excited on the flow conduit 1 . FIGURE 2 shows the flow conduit 1 in its undeflected position. FIGURE 3 shows 
a first elliptical radial mode where the term "first* describes the progression of complexity of the mode as it pertains to 

55 the number of waves above or below the undeflected position, along the length of the flow conduit 1 . For these exam- 
ples, it is assumed that the cross-sectional shape of the deflected conduit is elliptical, however there is a similar pro- 
gression in complexity of the natural modes pertaining to the cross-section as well, from elliptical to triangular to rec- 
tangular, and so forth for geometries of higher complexity. Assuming the cross-sectional deflected shape to be elliptical, 
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FIGURE 4 therefore represents a second elliptical radial mode, FIGURE 5 represents a third elliptical radial mode and 
FIGURE 6 represents a fourth elliptical radial mode. 

Any of the depicted mode shapes of FIGURES 3 through 6 and others not shown can be selected as the desired 
driven mode of vibration to cause the requisite Coriolis forces. In general, the resultant shape of the Coriolis force 
5 distribution, and thus the flow-related deflected shape of the conduit, will be similar to the next higher mode shape 
above the chosen driven mode shape. For example, for the preferred embodiment, the first elliptical radial mode is 
chosen as the desired driven mode of vibration as shown in FIGURE 3. The combination of the fluid flow 21 and conduit 
motion in the shape of the flow conduit 1 in FIGURE 3 cause a Coriolis force distribution and resulting conduit wall 
deflection that is similar in shape to the second elliptical radial mode as shown in FIGURE 4. Similarly, if the third 
10 elliptical radial mode of FIGURE 5 is chosen to be the driven mode of vibration, the resulting Coriolis force distribution 
and resulting deflected conduit shape will be similar to that of the fourth elliptical radial mode as shown in FIGURE 6. 

For a given conduit geometry, the magnitude of the flow induced deflection, as seen by the motion sensors 10 
through 1 5, is dependent on the mass flow rate and the frequency response of the conduit at the location of the motion 
sensors. Unfortunately, this frequency response is itself dependent on other factors (response variables) including (a) 
is the state of stress in the conduit from, for example, fluid pressure, end loading and gravity, (b) the elastic modulus of 
the conduit material (primarily a function of the materials temperature), (c) the density of the fluid in the conduit and 
(d) fluid viscosity. In addition, this frequency response is also dependent upon what type of motion detection method 
is used, i.e. velocity (as in the preferred embodiment), acceleration, displacement or time delay. 

it is the control of frequency response through specific geometric design or subsequent compensation (through 
20 either the application of a specific mathematical algorithm or direct determination of the response variables or coeffi- 
cients), that achieves the desired goal of insensitivity to these response variables. This is a central goal of the present 
invention. 

Four systems are disclosed that incorporate the desired insensitivity properties but vary in the complexity of either 
the sensor (mechanical portion of the device) or the electronics (signal processing portion of the device), since there 

25 js a natural tradeoff between the two. 

The first system to be discussed involves the most complex sensor geometry but the simplest signal processing 
necessary to achieve the desired results. Referring to FIGURE 7, the characteristic response curves 31 through 34 
for the conduit geometry of FIGURE 1 are plotted in graphical format. These response curves are calculated using 
finite element analysis ("FEA") methods to determine the change in the response of a particular sensor design to 

30 changes in fluid pressure, density or conduit elasticity, among other factors. 

A horizontal axis 36 of FIGURE 7 is essentially the progression of motion detection methods used to sense the 
Coriolis-induced conduit deflections starting with a time delay measurement 40 at the left, going through a displacement 
measurement 39, a velocity measurement 38 and an acceleration measurement 37 at the right end. Since the Coriolis- 
induced motion is a sinusoidal function of time, each successive measurement method along the horizontal axis 36 is 

35 mathematically related to the next by an integral (or differential) function of time and can be represented by integer 
multiples of the natural driven frequency of the conduit. For example, if accelerometers are used as the primary motion 
sensors on the conduit, integrating the acceleration -related output signals with respect to time or dividing the magnitude 
of these signals by the natural frequency Q of the conduit results in signals proportional in magnitude to the velocity 
of the conduit. Thus, velocity sensors instead of accelerometers can be used directly (as in the preferred embodiment). 

40 Similarly, integrating or dividing velocity signals by the natural frequency Q results in signals proportional to the dis- 
placement of the conduit. Thus, displacement sensors instead of velocity sensors can be used directly. Dividing dis- 
placement signals by the natural frequency Q results in signals proportional to the time delay of the conduit. Alternately 
time delay, rather than displacement, can be sensed directly. This has been the traditional method of choice. 

A vertical axis 35 of FIGURE 7 represents the change in the sensitivity of the device due to a given change in one 

45 of the response variables (i.e. pressure, temperature, density or stress). Accordingly, the pressure response curve 31 
shows that, with a given increase in fluid pressure (for example, increasing from 0 to 6.9 MNrrr 2 (1000 psi)), there is 
a positive increase in the sensitivity change of the device from left to right along the horizontal axis 36. Specifically, if 
time delay methods are used to measure the Coriolis-induced motion of the conduit, the working point along the hor- 
izontal axis 36 will be at a far left end 40 and, as pressure increases, according to the curve 31 the sensitivity of the 

so device will decrease (about -28% for the prescribed 6.9 MNnrr 2 (1000 psi) increase for this example). 

If displacement sensors are used to measure the motion, the working point along the horizontal axis 36 will be at 
the working point 39 where increases in pressure still cause a corresponding decrease in the sensitivity of the device, 
however not as much as if time delay methods had been used. Fortuitously, the pressure response curve 31 actually 
intersects a zero value at the working point 38 corresponding to the velocity of the conduit. Therefore, if velocity sensors 

55 are used as in the preferred embodiment, there will be no change in the sensitivity of the device with changes in fluid 
pressure. Finally, if accelerometers are used, the working point along the horizontal axis 36 will be at the working point 
37 resulting in an increase in sensitivity of the device with increases in fluid pressure. 

The density response curve 32 is similar to that just described for pressure, except it depicts the change in the 
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sensitivity of the device due to a given increase in the fluid density, for example increasing from the density of ambient 
air in the conduit to that of water. In general, the density response curves slopes in a direction opposite to that of the 
pressure response curves, thereby showing a decrease in sensitivity change as the working point moves from left to 
right along the horizontal axis 36. 

5 Again, if velocity sensors are used as the measurement method, there will be no change in the sensitivity of the 

device with changes in the fluid density. This fortuitous condition of having both the pressure and density response 
curves 31, 32 intersect a zero value at the same working point along the horizontal axis 36 and furthermore to have 
an intersecting working point 38 coincide with a fundamental measurement method such as velocity, is, of course, the 
desired result and is herein achieved by controlling the sensor geometric design parameters of flow conduit length, 

10 wall thickness, elastic modulus and the magnitude and location of individual masses such as the magnets used for 
both sensing and driving the conduit motion. 

The control of the slope and the intercepts of the pressure and density response curves 31 , 32 due to changes in 
specific geometric design parameters is hereinafter discussed in detail to allow the designer control over these char- 
acteristic curves. 

*5 The temperature response curve 33 of FIGURE 7 is shown sloping down toward the right indicating that an increase 

in the temperature (and thus a decrease in the elastic modulus) has a corresponding effect on the sensitivity of the 
device depending on which type of measurement method is used. The temperature response curve 33 intersects a 
zero value at the working point 38, corresponding to the velocity measurement method. Accordingly, if velocity meas- 
urement methods are used, changes in the elastic modulus due to temperature, will have little or no effect on the 

20 sensitivity of the device. In general, the shape of the temperature response curve 33 is less affected than are the 
pressure curve 31 or the density curve 32 by varying geometric design parameters. Since temperature is easily meas- 
ured and compensated for using simple electronic methods, the control of the shape and intercept of the temperature 
response curve 33 is usually not necessary. 

The stress response curve 34 of FIGURE 7 depicts the response of the device with an increase in axial tensile 

25 stress in the flow conduit of, for example, 138 MNm* 2 (20,000 psi). The stress response curve 34 accordingly shows 
that increases in axial tensile stress in the conduit result in a negative change in the sensitivity of the device, regardless 
of the type of detection method used. In general, the stress response curve 34 will not intersect a zero value within 
reasonable limits of design geometry. However, since axial tensile stress is easily eliminated by incorporating a bellows 
or slip pint arrangement in the conduit or measured using a strain gage, control of the shape and intercept of the axial 

30 stress response curve is usually not necessary. 

Therefore, for this example, a geometry is thus created that, by using velocity sensors, is unaffected by changes 
in pressure, density and temperature, simply by virtue of its geometry and motion detection method. Axial stress effects 
can also be eliminated by encompassing a bellows or slip joint arrangement onto the flow conduit, thereby making the 
device insensitive to all response variables except mass flow rate. An example of a set of geometric design parameters 

3S that meet this criteria are as follows: 



Flow Conduit Diameter 


2.5 cm 


0") 


Flow Conduit Length 


30.5 cm 


(12") 


Flow Conduit Wall 


1.2 mm 


(.049") 


Flow Conduit Metal 




316 SS (Elastic Modulus = 29e6) 


Drive Masses (2) 


5.24 g 


(.185 oz) (FIGURE 1 , items 16 and 17) 


Sensor Masses (4) 


1.16g 


(.041 oz) (FIGURE 1 , items 10, 11 , 14 and 15) 


Sensor Mass Locations 




25% & 75% of Conduit Length, Diametrically opposed 



45 

To use the geometry of this first embodiment, the sensor of FIGURE 1 is used in combination with the circuit of 
FIGURE 9, as follows. The requisite driven motion of flow conduit 1 is excited by the application of electrical signals 
of the appropriate phase and amplitude from a circuit component 81 into drive coils 16, 17 that are appropriately 
connected so that the resulting force causes the desired radial vibration, such as the first elliptical radial mode of 

so FIGURE 3. For maximum efficiency, the current supplied to coils 16, 17 is in phase with the desired velocity of the 
conduit motion. Once the desired vibratory motion is achieved, it is maintained at a constant level via an automatic 
gain control ("AGC") drive 80. 

With the driven motion established, signals from the diametrically opposed magnet/coil pairs 10, 11 near the inlet 
of the flow conduit 1 are summed together in a differential operational amplifier ("op-amp") 73 to yield signals propor- 

55 tional to radial motions and to eliminate signals proportional to common (bending mode) motion. Similarly, signals from 
the magnet/coil pairs 14, 1 5 near the outlet of the flow conduit 1 are summed together in a differential op-amp 74 to 
yield signals proportional to radial motions and to eliminate signals proportional to common (bending mode) motion. 
These resulting signals are then applied to both a summing op-amp 75 and a differential op-amp 76 to separate in- 
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phase and out-of-phase components. The output of the summing op-amp 75 is a sine wave proportional in amplitude 
to the driven velocity of the flow conduit 1 at the magnet/coil pair 10, 11, 14, 15 locations. Similarly, the output of the 
differential op-amp 76 is a sine wave proportional in amplitude to the velocity of the induced Coriolis deflections at 
those same locations. Since the Coriolis-induced velocity will normally be a small percentage of the driven velocity, 

s and since the in-phase portion of the signals from the differential op-amps 73, 74 may not match precisely in amplitude, 
the resulting difference from the differential op-amp 76 may contain a remnant amount of the in-phase component. 
This remnant value may be neglected if it remains constant over variations in operating conditions, or it may be sub- 
tracted by summing a countervailing, equivalent negative amount of the in-phase component through a component 
84. A demodulator 77 therefore receives input from the component 84 and determines the root mean square ("RMS") 

10 voltage or average voltage value of this sine wave using common techniques such as rectification, peak detection, 
asynchronous or synchronous demodulation. To maximize noise immunity, a synchronous demodulation technique 
using the driven frequency as the demodulation reference is the preferred technique. If the strict geometric design 
parameters as just described are adhered to and if axial stress is held to be negligible, no further compensation is 
required for this signal, since as previously described, it will not vary with pressure, density or temperature variations. 

is The demodulator 77 therefore can provide an output signal (M1) 82 that is proportional to the mass flow rate and 
insensitive to these other variables. This represents the first embodiment that imposes the most stringent requirements 
on geometric design but has the simplest signal processing circuitry. 

Because of normal material and manufacturing variations, the response curves 31 through 33 may not intersect 
perfectly as shown on FIGURE 7. As will be discussed later in detail, these response curves can be controlled by 

20 design. However, after a sensor has been manufactured and found to be slightly "off" due to these variations, a final 
correction can be made during a subsequent calibration procedure. To effect this correction, the sensor is flow-tested 
at two different pressures and densities to determine the slope and intercept of the response curves. With this infor- 
mation, additional masses can easily be added at appropriate points along the length of the flow conduit 1 to shift the 
response curves back to the correct intercepts. 

25 The geometric design requirements of this first embodiment may impose restrictions upon the designer that are 

undesirable. As a practical matter, there are many possible sensor attributes that a designer may want to optimize for 
a particular design. Among these are pressure drop with flow, sensitivity, natural frequency, safety and corrosion re- 
sistance. Each of these can be effected by geometric design parameters and it is therefore desirable to allow the 
designer as much freedom as possible to optimize the most important attribute for a given design. 

30 The second system therefore addresses this problem and allows the designer more freedom to optimize geometric 

design parameters such as flow conduit length, wall thickness and magnet size and location, without the restriction of 
having the pressure and density response curves 31 , 32, intersect at a primary measurement method such as velocity 
38. 

As an example of a sensor geometry for the second system the following design parameters are representative: 

35 



Flow Conduit Diameter 


2.5 cm 


0") 


Flow Conduit Length 


25 cm 


(10") 


Flow Conduit Wall 


0.71 mm 


(.028") | 


Flow Conduit Metal 




316 SS (Elastic Modulus = 29e6) ! 


Drive Masses (2) 


2.3 g 


(.08 oz.) each (FIGURE 1, items 16 and 17) 


Sensor Masses (6) 


1-16 g 


(.041 oz) each (FIGURE 1, items 10 through 15) 


Sensor Mass Locations 




25%, 50% & 75% of Conduit Length, Diametrically opposed 



45 This second embodiment is therefore shorter, has thinner walls and lighter drive masses than the first embodiment, 

which may be desirable to optimize some of the aforementioned sensor attributes such as pressure loss and sensitivity. 

Referring to FIGURE 8, pressure and density response curves 51 , 52 are representative of the response charac- 
teristics of this second system. These curves indicate similar shape and slopes to those of the first embodiment (the 
curves 31 , 32 of FIGURE 7), however, the significant difference is that their zero value intercepts 55 occur at a working 

so point along an horizontal axis 56 that does not coincide with a primary measurement method. However, due to the 
mathematical relationship between the measurement methods as earlier described, it is not necessary to achieve 
coincidence with a measurement method. No matter the type of measurement method chosen for the primary motion 
sensing, the resulting signals can be mathematically transformed to any arbitrary working point along the horizontal 
axis 56 by applying to the resulting flow-related signals, a function proportional to Q n , where Q is the natural frequency 

55 of the conduit at the time the flow signals are attained and n is the exponential number (either positive or negative) 
necessary to achieve the desired transformation. In this example, velocity signals were the chosen measurement 
method using a common magnet and coil arrangement to derive signals proportional to the velocity of the Coriolis- 
induced motion. 
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Dividing the magnitude of these velocity signals by Q 1 equates to using displacement sensors and reduces the 
pressure and density sensitivity by the amount shown on the pressure and density response curves 51 , 52, however 
this does not eliminate the problem. Dividing these velocity signals by O 2 equates to using time delay measurement 
and, for this example, reduces the pressure and density effects (and reverses their polarities) but, again, does not 
5 eliminate them. 

The use of integer real values (either positive or negative) for the exponential number n functionally transforms 
the fundamental sensor signals from one fundamental type (velocity) to another (displacement, etc.) . While commer- 
cially-available transducers are usually restricted to these fundamental types corresponding to integer real values of 
n, the exponent number in the algorithm is not restricted to integer real values and, indeed, the correct function for this 

io geometry is to divide the magnitude of the flow-related velocity signals by n 1 6 (a noninteger number for n). This renders 
the resulting flow signals insensitive to both pressure and density at the same time. Physically, this equates to using 
a hypothetical sensor providing signals proportional to a blend of 60% time delay and 40% displacement. This can 
actually be achieved by using these two methods and blending their signals in the correct proportions or by using 
displacement sensors and applying a partial integration of the signal. However, these methods are analogous and tend 

is to add unnecessary complexity to the simple method described herein for the preferred embodiment. 

A temperature response curve 58 and an axial stress response curve 59 are similar to the corresponding temper- 
ature and axial stress response curves 33, 34 of FIGURE 7 of the first embodiment. These curves 58, 59 were effectively 
unchanged by the geometric difference between the first and second embodiments. By dividing the flow-induced ve- 
locity signals (M1 from the demodulator 77 of FIGURE 9) by Ci 1 - 6 , the effective working point along axis 56 of FIGURE 

20 8 is transferred from velocity to an intercept working point 55 that eliminates the effects of both pressure and density. 
However, unlike the first embodiment, it creates a sensitivity to temperature, as can be seen by the positive value of 
the temperature response curve 58 above the intercept working point 55. Therefore, working at the working point 55, 
the sensitivity of the device increases with an increase in temperature. Similarly, the axial stress response curve 59 
indicates that, at the working point 55, there will be a negative effect on the sensitivity of the device with an increase 

25 in axial stress. 

To use the algorithm and compensate for both temperature and axial stress, the circuit of FIGURE 9 incorporates 
all the components used for the first embodiment and further incorporates additional circuit components that were not 
used for the first embodiment. 

To compensate for the effect of temperature, a temperature sensor 18 and a corresponding scaling component 70 

30 are used to create a scaled compensation value appropriate for the chosen working point 55 along the horizontal axis 
56 in combination with the temperature response curve 58. Similarly, if axial strain is not held to be negligible or constant, 
its value can be measured and compensated for via the strain gage 19 and a corresponding scaling component 71, 
used to create a scaled compensation value appropriate for the chosen working point 55 along the horizontal axis 56 
in combination with the axial stress response curve 59. These compensation values from the scaling components 70, 

35 71 , can be summed and further scaled in the compensation component 72 and then applied to the computational 
component 79. 

The computational component 79 receives an uncompensated mass flow rate signal (M1) from the demodulator 
77, stress and temperature compensation values B from the compensation component 72 and the value of the driven 
frequency £2 from the component 78. These values are then combined within the computational component 79 to apply 
40 the following algorithm: 

M2 = M1 *B*(A/fl°) (1) 

45 Where: 

M2 = Final compensated flow rate signal; 

M1 = Uncompensated flow rate signal; 

B = Temperature and strain compensation values; 

50 A = Scaling, i.e. [0-pressure frequency(Hz)] / [0-pressure sensitivity(Volts/Flow Rate)]; 

CI = natural frequency; and 

n = Exponent appropriate for the working point along axis 56 (1 .6 for this example). 

The resultant output signal (M2) 83 from the computational component 79 is therefore insensitive to pressure and 
55 density changes by virtue of the (ft n ) portion of the algorithm and is further insensitive to temperature and axial stress 
changes by virtue of the (B) portion of the algorithm. As a practical matter, the functionality of the computational com- 
ponent 79 can be implemented via analog components such as Burr Brown part number 4302, or Analog Devices part 
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10 



15 



number AD538. Similarly, this functionality can easily be converted to a digital design with the algorithm being imple- 
mented as a sequence of computer instructional codes executable in a general purpose computer in a conventional 
manner. 

FIGURE 10 depicts uncompensated actual data from flow testing of a sensor designed according to this second 
embodiment. The horizontal axis represents the flow rate in pounds of air per hour and the vertical axis represents the 
percentage deviation from the reference meters. Data sets 102 through 105 represent the percentage deviation of the 
device from the reference meters, at various flow rates and at 4 different pressures. The data set 102 was taken with 
air flow at 100 psi pressure. Similarly, the data sets 103, 104, 105 correspond to pressures of 1 .4, 2. 1 and 2.8 MNnr 2 
(200, 300 and 400 psi), respectively. These uncompensated data sets correspond to signal (M1) 82 of the circuit of 
FIGURE 9. This corresponds to working at the velocity working point along the horizontal axis 56 of FIGURE 8 that 
predicts the positive effect on sensitivity with increasing pressure as shown by the positive value of the pressure re- 
sponse curve 51 above the velocity working point. Therefore, as the pressure was increased for each data set, the 
meter sensitivity increased causing increasingly higher deviation values as seen in FIGURE 10. 

The temperature and axial stress values were held to be constant during the tests, eliminating a need to include 
compensation for these parameters in the algorithm. The following Table I represents the observed operating frequency 
of the tested device as a function of pressure. These frequency values were then incorporated into the algorithm with 
the results shown in FIGURE 11. 



TABLE I 



20 



25 



psi 


Frequency (Hz) 


0 


2951.2 


100 


3014.7 


200 


3075.2 


300 


-3135.7 


400 


3196.3 



30 



35 



40 



45 



SO 



55 



Data sets 110 through 113 of FIGURE 11 therefore depict the percentage deviation with respect to the reference 
meters after the application of the algorithm and corresponds to the compensated output signal (M2) 83 of FIGURE 
9. FIGURE 1 1 shows that the application of the algorithm, by virtue of the increasing frequency term in the denominator 
with increasing pressure, effectively offsets the increasing sensitivity of the velocity signals (M1 ) 82 by the exact amount 
necessary to render the final output signals (M2) 83 insensitive to both pressure and density changes. 

This second system therefore constitutes an system that is inherently insensitive to both pressure, density, tem- 
perature and axial stress by virtue of its geometry in combination with a compensation algorithm. This embodiment 
has the tremendous advantage of allowing the designer a wide range of geometric freedom to maximize any desired 
performance attribute. However, this embodiment requires more sophisticated signal processing than does the first 
embodiment due to the application of the requisite algorithm. 

A third system allows more geometric design freedom than the previous two embodiments without adding more 
complexity to the electronics. The tradeoff to this embodiment is that it is insensitive to pressure and density changes, 
but only for specific fluid properties, specifically fluids with a predictable relationship between their pressure and density. 
FIGURE 12 depicts pressure and density response curves 151, 152, respectively, according to this third embodiment. 
By allowing the designer complete geometric freedom, sensor geometries can be created having pressure and density 
response curves 151, 152 that do not intersect zero values (155, 156) at the same working point along a horizontal 
axis 157. This situation allows the designer maximum geometric design freedom but presents the problem of how to 
compensate for changing response variables. As a practical matter, zero intercepts 155, 156 can generally be main- 
tained acceptably close to one another for most sensor geometries by design. However, both geometric design pa- 
rameters, as well as material and manufacturing anomalies, can alter the location of these intercepts and thereby cause 
them to be non coincident. 

One solution to this problem according to the third system is specific to the intended fluid to be measured. For 
example, if it is known that the fluid to be measured will not appreciably change density over the intended operating 
conditions, then the working point is chosen at the zero intercept 155 along the axis 157 and the appropriate value for 
the exponential number n is applied using the aforementioned algorithm. For this example, to move the working point 
from velocity to the working point zero intercept 1 55 requires a value for n of approximately 1 .3. Similarly, if the appli- 
cation is known to be at constant pressure, then a working point should be chosen at the zero intercept 156 and the 
appropriate value for the exponential number n should be applied using the aforementioned algorithm. For this example, 
to move the working point from velocity to the zero intercept 1 56 requires an exponent value for n of approximately 1 .7. 
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If there is a known relationship between the pressure and the density of the fluid to be measured, for example, by 
using the ideal gas equation for a gaseous fluid, then a third, net resultant response curve 1 53 can be plotted in FIGURE 
1 2 that is the net resultant response curve for the combined effect of both pressure and density for that particular fluid. 
To compensate for this net effective response curve 153, a working point 154 is used and the appropriate value for the 

5 exponential number n should be applied using the algorithm of equation (1). For this example, to move the working 
point along the axis 157 from velocity to working point 154 requires a value of approximately .6 for the exponential 
number n in equation (1). Since there may be a temperature dependence between the pressure and density of a given 
fluid, the temperature compensation value (not shown) can be modified in the scaling component 70 and compensation 
component 72 of FIGURE 9 to accommodate this dependence. 

io The slope of the pressure response curve 151 is more affected by changes in the wall thickness of the conduit 

than is the density response curve 152. Therefore, by adjusting the flow conduit wall thickness for a particular design 
(in this example, an increase in wall thickness), the slope of the net resultant response curve 153 can be made to be 
nearly coincident with and parallel to the horizontal axis 1 57, thereby allowing the designer to work virtually anywhere 
along the horizontal axis 157 without errors resulting from either pressure or density for that particular fluid. As a 

is practical matter, the net resultant response curve 153 is not, in general, linear and is therefore not coincident with the 
horizontal axis 157 over its entire length. However, the resulting error may be acceptable for many practical designs 
and will allow the designer to pick a convenient value for the exponential number n (either positive or negative, integer 
or noninteger) in the algorithm of equation (1 ) or an integer real value to allow the use of a fundamental measurement 
method (displacement etc.). 

20 a fourth system involves the determination of pressure, density or both and subsequent compensation for either 

or both of their effects on any arbitrary flow conduit geometry. The preferred method to determine the pressure, density 
or both for this embodiment is to vibrate the sensor in two natural modes of vibration simultaneously, with the restriction 
that the vibrational response of the device to changes in either pressure or density must be different for the two different 
modes. This difference allows for simultaneous solution of the pressure and density and subsequent compensation 

25 according to the device's characteristic curves. This embodiment has no restrictions on where the pressure or density 
sensitivity curves intersect zero or each other, thereby allowing maximum freedom to the designer to optimize any 
chosen design attribute such as conduit length, conduit walls thickness or magnet weight and to apply that geometry 
to any fluid to be measured. 

FIGURE 1 3 depicts pressure and density sensitivity curves 161 , 162 respectively for a representative flow conduit 
30 geometry. FIGURE 13 depicts the pressure and density response curves 161 and 162 having zero intercepts 165, 167 
that are not coincident and an intersection 168 that is not a zero value. As previously mentioned, the zero intercepts 
165, 167 are normally close to each other, however it is possible to get large separations between the two by design. 
Therefore FIGURE 13 serves as a good example for the fourth embodiment where the slopes and intercepts of these 
response curves are not restricted. 
35 Using the simultaneous solution method (described in detail hereinafter) the value of pressure and density are 

both determined and their effects compensated for. For example, referring to FIGURE 1 3, if velocity signals are used 
for motion detection, then the necessary value for pressure compensation is the negative of the working point 163 of 
the curve 161, or approximately -30% per 1000 psi (since the response is positive with pressure, the compensation 
should be negative). Similarly the working point 164 on the curve 162 indicates a density compensation value of 8% 
40 per 1 gram/cc. Using the simultaneous solution method, the values of pressure and density are both accurately deter- 
mined and then applied in conjunction with the aforementioned compensation values to correct the resultant signals 
for both pressure and density. For this example, if the pressure was determined to be 1 .7 MNm 2 (250 psi), then the 
pressure compensation value is: 

45 -2 -2 

1 .7 MNm (250psi) * (-30%/6.9 MNm (1000psi)) = -7.5% compensation value. 

Similarly, if the density was therefore determined to be .5 gm/cc then the density compensation value is: 

so 

.5g/cc * (8%/1g/cc) = 4% compensation value. 

To use the simultaneous solution method, for determination of both pressure and density, two modes of vibration 
are excited on the sensor at the same time, with the requirement that their natural frequencies must change as a 
ss function of pressure or density at different rates. 

FIGURE 14 is an isometric view of a three axis graph depicting a pressure (right axis), density (left axis) and 
frequency (vertical axis), relationship for two radial modes of vibration for the sensor geometry of FIGURE 1 . A surface 
180 depicts frequency values for various combinations of pressure and density for the first elliptical radial mode of 
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vibration. Similarly, a surface 181 depicts frequency values for various combinations of pressure and density for the 
third elliptical radial mode of vibration. The surfaces 180, 181 are mathematical representations and are determined 
in a compensation factor determination process by subjecting a vibrating sensor to various combinations of pressure 
and density to determine a representative number of frequency values. Visual inspection of the surfaces 180, 181 

s indicate that, for both modes, frequency increases as pressure increases and decreases as density increases, but in 
different amounts. Other factors that affect these natural frequencies include temperature, conduit stress and viscosity. 
However, these other variables can be separately measured via temperature sensors, strain gages, drive power meas- 
urements or other similar means and are therefore not included as independent variables in the simultaneous solution. 
During operation, both modes of vibration are excited and their frequency values determined. Knowing the fre- 

10 quency of each mode then allows a mathematical determination of a plane of constant frequency intersecting each of 
the surfaces 180, 181. The intersection of a plane of constant frequency and either of the surfaces 180, 181 thereby 
defines a curve of constant frequency 182 thereon. The constant frequency curve 182 is therefore representative of 
all the possible combinations (within the limits of the graph) of pressure and density that causes that particular frequency 
for the first elliptical radial mode of vibration. Similarly, a constant frequency curve 183 is representative of all possible 

is combinations of pressure and density that cause that particular frequency for the third elliptical radial mode of vibration. 

FIGURE 15 is another view of the graph of FIGURE 14, but viewed parallel to the frequency (vertical) axis and 
perpendicular to the pressure/density plane. Because of the different shapes of the surfaces 1 80, 181 of FIGURE 14, 
the curves 182, 183 are not parallel to each other and, in this view, intersect each other at a solution point 184. The 
solution point 184, therefore represents the only combination of pressure 185 and density 186 that can cause the 

20 aforementioned two frequencies to occur at the same time. Therefore, by solving for this intersection working point, 
the pressure 185 and density 186 can be simultaneously determined and their values used for appropriate compen- 
sation as previously mentioned. 

This system is therefore rendered insensitive to pressure and density by virtue of its geometry in combination with 
direct solution of both pressure and density by the method of simultaneous solution and direct application of appropriate 

25 compensation factors. The system therefore represents the highest degree of geometric design freedom to optimize 
any desired attribute, but requires the most sophisticated signal processing to effect the simultaneous solution. 

FIGURES 16 and 17, taken together, depict a signal processing circuit that can be used in conjunction with the 
sensor geometry of FIGURE 1 for this fourth embodiment. FIGURES 16 and 17 include many of the circuit components 
used previously in FIGURE 9 and further include additional components related to sensing and driving the requisite 

30 secondary mode of vibration and effecting the simultaneous solution. 

To use this fourth system, the geometry of FIGURE 1 is excited in two modes of vibration simultaneously, the first 
mode preferably being the first elliptical radial mode (as shown in FIGURE 3) and the second being the third elliptical 
radial mode (as shown in FIGURE 5). These vibrations are sensed via the magnet/coil pairs 10 through 15. However, 
for this embodiment, the magnet/coil pairs 12, 13, being located midway down the length of flow conduit 1 , sense the 

35 motion of the third elliptical radial mode (as shown in FIGURE 5) in a phase opposite (180°) from that sensed by the 
other magnet/coil pairs 10, 11, 14, 15. This 180° phase inversion is used by summing and differential op-amps 201, 
202 of FIGURE 1 6 to sum appropriate amounts of the signals received from each magnet/coil pair to separate the first 
elliptical radial mode Q from the third elliptical radial mode a. The summing op-amp 201 sums signals from all 6 magnet/ 
coil pairs, reinforcing the £2 signals, since they are all in the same phase. However, for the second mode ft, signals 

40 from the magnet/coil pairs 12, 13 are 180° out of phase compared to those of 10, 11, 14, 15. Therefore, by amplifying 
the signals appropriately from the differential op-amps 73, 74, 200, a cancellation of O-related signals is accomplished 
in the summing op-amp 201 , thereby isolating the ^-related signals. Similarly, using the same technique, the differential 
op-amp 202 thereby reinforces Q-related signals and cancels Q-related signals. By isolating these two signals from 
each other, each can be separately controlled for amplitude purposes and each can be accurately measured for use 

45 in obtaining a simultaneous solution. Q and a signals from the components 201 , 202, respectively, are both delivered 
to an AGC drive 80 and a power driver 81 to reinforce and control each mode's amplitude. In addition, Q and Q signals 
from the summing and differential op-amps 201, 202 are sent to analog to digital ("A/D") frequency converters 206, 
207, respectively, that convert the signals from analog sine waves to digital numbers representing their frequencies. 
These numbers are then sent to a microprocessor 220 that uses these values to develop the simultaneous solution. 

50 To use the simultaneous solution method, the pressure-versus-density-versus-frequency characteristics of a par- 

ticular design must be determined to define mathematically the equations for the surfaces 180, 181. This is most 
accurately done through a compensation factor determination procedure by subjecting the sensor to various pressures 
and densities and noting the resulting frequencies for both Q and Q. With a number of specific frequency working points 
known, both surfaces 180, 181 can be mathematically defined as a function of frequency as follows. 

55 

n = F1 Q + (C1)*P - (C2)*D - (C5)T + (C6)e - (C9)v (2) 
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a = F2 0 + (C3)*P - <C4)*D- (C7)T + (C8)e - (C10)i> (3) 



20 



25 



30 



35 



40 



45 



50 



55 



Where: 



ft.Q Measured Frequency for first and second modes of vibration, respectively (Hz); 

F1 0 , F2 0 Zero pressure/density frequency values for both modes, respectively (Hz); 

C1 - C10 Constants determined for each sensor or design. (Hz/psi, Hz/g/cc etc.); 

P Fluid Pressure (psi); 

10 D Fluid Density (g/cc); and 

T.e.-o Temperature, Axial stress, Viscosity. These variables can be separately measured via the sensor 18 and 
strain gage 19 and need not be carried as independent variables in the simultaneous solution. 

For the given sensor geometry for this fourth system, the actual values for the pressure and density-related coef- 
15 ficients are similar to the following: (The temperature, stress and viscosity related terms can be separately measured 
and added as constants to the following equations, but for simplicity are not shown here.) 



a = 2951 + (.7)*P - (16.65)*D (4) 

ft = 4408+ (.51 )*P - (25.9)*D (5) 

These two equations therefore define the three dimensional surfaces 180, 181 in terms of pressure, density and 
frequency for both modes of vibration. It is important to note that the coefficients preceding the pressure and density 
terms of equations (4) and (5) above have different values thereby according different slopes to the surfaces 1 80, 1 81 . 
These differences are necessary to cause the curves 182, 183 to intersect in the plane of FIGURE 15 and thereby 
allow a simultaneous solution. The larger the slope differences between surfaces 180, 181, the more divergent the 
curves 182, 183 are and the more accurately their solution point 184 can be determined. 

Next, the meter is placed in operation and the actual in-situ Q. and Q frequencies are determined in the A/D fre- 
quency converters 206, 207 and their values delivered to the microprocessor 220. This reduces equations (4) and (5) 
above to functions of pressure and density only and in effect define the two curves 182, 183 of FIGURES 14 and 15. 

Next, the one unique solution to these two equations (pressure 185 and density 186) are determined in the micro- 
processor 220 using the Newton-Raphson method of simultaneous solutions. This method uses the two equations (4) 
and (5) and their derivatives to solve for two equations with two unknowns where the equations can be non-linear. 
Since the equations (4) and (5) above are in polynomial form, their derivatives are easily determined. While equations 
(4) and (5) above are linear in nature, the actual surfaces 180, 181 will have a nonlinear component. For higher accu- 
racies, these nonlinear components can also be added to equations (4) and (5) since the Newton-Raphson method 
can solve for nonlinear equations as well. In lieu of supplying the equations and derivatives in polynomial form, tables 
of values can be supplied in the form of a look up table. This method simplifies necessary curve fitting to reduce the 
data describing the surfaces 180, 181 into polynomial form and lends itself to accurately describing the nonlinear 
characteristics of the surfaces 180, 181 . 

For a reference to the described simultaneous solution method and a sample of computer code to implement it, 
refer to the book "Numerical Recipes in C", Cambridge University Press, second edition, ISBN 0-521-43108-5, § 9.6 
entitled, "Newton-Raphson Method for Nonlinear Systems of Equations." Using this method in conjunction with equa- 
tions (4) and (5), the solution point 184 can be determined and its corresponding values for pressure 185 and density 
186 are thereby determined. 

The circuit of FIGURE 1 6 further uses temperature signals from the scaling component 70 and strain-related signals 
from the scaling component 71 , obtaining the signals via a multiplexer 203 and an A/D converter 204. In addition, the 
multiplexer 203 also receives signals from the differential and summing op-amps 76, 201 , 202 to allow the microproc- 
essor 220 to control various other aspects of the circuit, for example, the gain settings of the scaling and compensation 
components 70, 71 , 73, 74, 200, as well as controlling the amplitudes of both modes of vibration. Flow-related velocity 
signals from the differential op-amp 76 are converted to digital values 209 in the A/D converter 205 and sent to the 
microprocessor 220. Therefore, the microprocessor 220 receives the values of 0, 21 0 and a 21 1 from the A/D frequency 
converters 206 and 207, the value of flow-related velocity signals 209 from the A/D converter 205 and the values of 
temperature, strain, O. and Q driven velocities and others 208 from the A/D converter 204. In addition, the microproc- 
essor 220 receives coefficient data defining the characteristics of the sensory (coefficients 1 through 10 of equations 
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(2) and (3) above) from the storage component 222. These signals therefore comprise the necessary information to 
effect the simultaneous solution and render the necessary compensation. 

Referring now to FIGURE 1 3 f with the pressure and density values now determined, the correct compensation 
values corresponding to the working points 163 (for pressure) and 164 (for density), can be directly applied in the 
s microprocessor 220. In addition, the microprocessor 220 compensates for temperature, axial stress and viscosity, as 
necessary. 

Finally, the microprocessor 220 supplies information to an output module 221 that can create output signals to the 
user including mass flow rate, pressure, density, temperature, viscosity and other calculated values. These output 
signals can further be used to control process equipment coupled to the mass flow meter, such as pumps to maintain 
10 a proper flow rate; valves to maintain a proper fluid density or viscosity; or temperature regulators to maintain a proper 
temperature. 

Under certain circumstances, the simultaneous solution method may prove more accurate for solving for pressure 
than density or vice versa. An example of this is highly incompressible gasses, such as hydrogen, where the density 
is very low and changes little with large changes in pressure. Therefore, as a modification to this fourth embodiment, 

is higher accuracies may be attained by a combination of applying the previously-defined algorithm of equation (1 ) of the 
second embodiment to eliminate the effects of one parameter, for example density, and then directly compensating for 
the other parameter, for example pressure, using the pressure value determined by the simultaneous solutions. This 
greatly reduces the accuracy requirement of the simultaneous solution since no compensation for density is needed 
and a much smaller compensation for pressure is used. 

20 For this example, where the pressure and density response curves are as shown in FIGURE 13, the flow-related 

velocity signal (M1) 209 of FIGURE 16 is divided by a function proportional to (Q) 7 to shift the working point from 
velocity to a working point 165, to eliminate the effect of density alone. Then, using the determined value of pressure 
from the simultaneous solution, a pressure-compensation value is directly determined from a working point 166 on the 
curve 161 and the appropriate compensation value directly applied, or approximately 20% per 1 000 psi for this example. 

25 This fourth system therefore is inherently insensitive to pressure, density, temperature and axial stress by virtue 

of its geometry in combination with the algorithm of equation (1 ) and direct compensation using the determined values 
for pressure, density or both from the simultaneous solution and direct temperature and conduit strain measurements. 

The four systems just described eliminate the effects of pressure, density, temperature and axial stress by dealing 
with the response characteristics of a given sensor design with different signal processing methods. In the first system, 

30 the geometry was designed to have the pressure, density and temperature response curves intersect each other and 
zero coincident with a fundamental measurement method (velocity). The second system forces the pressure and den- 
sity response curves to intersect zero at the same working point but not coincident with a fundamental measurement 
method, etc. For a designer to decide which system is best for a particular application, the designer needs to know 
how to control the slope and zero intercept of the pressure and density response curves for a prospective design. The 

35 following discussion therefore gives guidelines with respect to controlling these response characteristics. 

The slopes and zero intercepts of the response curves in FIGURES 7, 8, 12 and 13 are manipulated by adjusting 
the mass and stiffness properties of the vibrating portion of the conduit assembly. For the sensor geometry of FIGURE 
1 , the manipulation of stiffness properties is accomplished by modifying the conduit diameter, wall thickness, elastic 
modulus and conduit length, while the manipulation of mass properties is accomplished by modifying the conduit density 

40 and the number, size and location of individual lumps of mass on the conduit. 

While individual lumps of mass can be non-functioning, the necessary mass at a location can usually be combined 
with the functionality of a motion driver or motion detector by making the mass a magnet, a coil or some other type of 
motion driver or detector. It is important to realize that the effects of manipulating these design parameters are highly 
interrelated. For example, increasing the wall thickness of the conduit not only increases the stiffness to radial vibrations, 

45 jt also increases the mass of the conduit uniformly and diminishes the effect of individual lumps of mass located on 
the conduit. Therefore, keeping in mind that changing one aspect of the sensor geometry effects many others and that 
the effect of a single change is not absolutely deterministic without analysis, general trends can be identified and will 
now be discussed giving the designer sufficient capability to control response curve shapes and intercepts. 

In general, a single straight flow conduit without any extra masses placed upon it has pressure and density response 

50 characteristics similar to those depicted in FIGURE 18. This figure compares the pressure and density response char- 
acteristics of three 2.5 cm (1 ") diameter, 0.71 mm (.028") wall thickness stainless steel flow conduit geometries without 
any additional lumped masses. Data sets 230, 233 are the pressure and density response curves for an 20 cm (8") 
long conduit, data sets 231 and 234 are for a 25 cm (10") long conduit and data sets 232, 235 are for a 30 cm (12") 
long conduit. One important aspect of the data sets 230, 231, 232, 233, 234, 235 of FIGURE 18 is that their zero 

55 intercepts all closely converge on the time delay signal processing method on the left. This is a natural characteristic 
of a perfectly round straight conduit with uniform mass distribution (no additional masses) and uniform stiffness distri- 
bution (no added stiffness e.g. from brazed-on magnet hardware or changes in the wall thickness). If a bare conduit 
can be excited into the requisite mode of vibration and that motion sensed without adding any localized masses or 
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stiffness (such as magnets) then the working point on the horizontal axis of FIGURE 18 is the time delay method at 
the far left. As a practical matter however, to cause the requisite vibration and sense its motion, motion drivers and 
sensors are normally employed that add localized mass and stiffness properties to the conduit that alters these response 
curves away from converging on time delay. 

s To exemplify the effect of adding two 5.24 g (.185 ounce) masses at the center of a conduit similar to the magnet/ 

coil pairs 16, 17 on the flow conduit 1 of FIGURE 1, the resulting response curves for the three conduit lengths of 
FIGURE 1 8 are plotted in FIGURE 1 9. Comparing the bare conduit response curves of FIGURE 18 to those in FIGURE 
1 9 indicates that the addition of center mass reduces the slopes of the response curves and shifts their zero intercepts 
to the right of the time delay method by varying amounts depending on the conduit length. 

10 FIGURE 20 depicts the response characteristics of the same three conduit lengths and center masses of FIGURES 

1 9 but with an increase in wall thickness to 1 .2 mm (.049"). The effect of this increase in wall thickness further reduces 
the slopes of the response curves and further shifts the zero intercepts of the response curves to the right. 

The effect of adding outer masses; such as the magnet/coil pairs 10, 11 , 14, 15 on the flow conduit 1 of FIGURE 
1 , varies depending on how much the pairs 10, 11 , 14, 15 weigh in comparison to center masses of the drive coils 16, 

is 17 and their position along the length of flow conduit 1. If the outer masses are located close to the center masses, 
the outer masses exacerbate the effect of the center masses. Alternately, if the outer masses are located near the 
conduit ends (near the manifolds 4, 5 of FIGURE 1) the effect of the outer masses is negligible. If the outer masses 
are located midway between the conduit ends and the center masses, the outer masses detract from the effect of 
center masses. 

20 Therefore, the following Table 1 1 summarizes the general effects of increasing any one particular design parameter 

as it pertains to the slope of the response curve, its zero intercept, the flow conduit frequency and the sensitivity. 



TABLE II 



25 



Design Parameter 
(Increases) 


Response Curve 
Slope 


Response Curve 
Intercept 


Flow Conduit 
Frequency 


Sensor Sensitivity 


Conduit Length 


Negligible 


Shifts Right 


Decreases 


Increases 


Wall Thickness 


Decreases 


Varies Negligibly 


Increases 


Increases 


Elastic Modulus 


Decreases 


Varies Negligibly 


Increases 


Increases 


Center Mass 


Decreases 


Shifts Right 


Decreases 


Decreases 


Outer Mass 


varies 


varies 


Decreases 


Decreases 



3S As previously mentioned, the effects of changing any one design parameter are interrelated with the effects of the 

other parameters and, therefore, the final result of a single change is not absolutely deterministic without analysis. 

The methods and algorithms delineated in the specification herein have been applied to Coriolis flow sensors with 
a single conduit vibrating in a radial mode of vibration, since that is a preferred embodiment. However, the control of 
the slope and intercepts of the various response curves and the algorithm of equation (1 ) also pertain to traditional 

40 bending mode type Coriolis flow meters as well. It is outside the scope of the Detailed Description to expound upon 
the effects of changing design parameters as they pertain to traditional bending mode Coriolis type flow sensors since 
there are currently many complex conduit shapes and geometries available ranging from simple straight conduits 
vibrating in the first bending mode to Q-shaped conduits vibrating in second order bending or twist modes. 

Given all the possible design parameters for bending mode conduit geometries, it is impractical to construct a 

4S simplified summary of their effects on response as in Table II, above. However, the basic principles are the same in 
that by appropriate application of localized mass and stiffness properties upon the flow conduit, the response curves 
of a given bending mode sensor can be controlled so that the zero intercepts coincide with each other at some working 
point along the horizontal axis, wherein the appropriate value of n for the exponent of equation (1 ) can be applied to 
correct for pressure and density at the same time. Alternately, if the response curves do not intersect at a zero value, 

so the method of simultaneous equations as explained above for radial mode meters can be applied to bending mode 
meters to calculate directly the pressure and density and directly compensate for those errors. 

As an example, FIGURE 21 depicts pressure and density response curves 260, 261 for a 10" long, .049" wall 
thickness stainless steel conduit with no masses on the conduit and operating in the first bending mode of vibration. 
Since there are no added localized masses or stiffnesses on this example, it is comparable to the graph of FIGURE 

ss 18 (a radial mode meter with no added masses). 

Comparing FIGURES 21 to 18 indicates that, in general, a bare flow conduit without added masses or stiffnesses 
has response curves that tend to converge near the time delay method on the left side of the graphs. Another obser- 
vation that can be made between FIGURES 18 and 21 is that the slope of the pressure response curve 260 in the 
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bending mode is significantly less than the slope of the curves 230 through 232 operating in the radial mode. This 
slope difference is due to the fact that bending mode vibrations are not as greatly affected by internal pressure as are 
radial modes of vibration. However, bending mode conduit designs having thinner walls and more convoluted shapes 
than this example, have increasingly higher slopes on their pressure response curves. Therefore, the error caused by 

5 pressure and density changes can be significant on bending mode meters and can be corrected using the methods 
described herein for radial mode meters. 

As with radial mode meters, localized masses and stiffnesses are normally added to a bare flow conduit operating 
in a bending mode of vibration as necessary for motion drivers and sensors. The effects of these added masses and 
stiffnesses are similar to those on radial mode meters. FIGURE 22 therefore is a final example of the applicability of 

to the present invention as it applies to bending mode meters. Pressure response curves 270 through 272 and density 
response curves 273 through 275 indicate the response of conduit lengths of 20 cm (8 a ), 25 cm (10") and 30 cm (12"), 
with 0.71 mm (.028") wall thickness, of stainless steel, operating in the first bending mode of vibration and with two. 
185 ounce drive masses at the center of the conduit length. These response curves show the susceptibility of bending 
mode meters to localized masses and stiffness and show that the algorithm of equation (1 ) can be advantageously 

is applied on bending mode meter designs as well. For the example of FIGURE 22, using the 20 cm (8") long conduit 
having the corresponding response curves 270, 273 and assuming velocity sensors were used to acquire the flow 
signals, an exponent value of approximately 1.8 in equation (1) can be used to correct for most of the pressure and 
density-related errors on these signals. 

Although the present invention and its advantages have been described in detail, those skilled in the art should 

20 understand that they can make various changes, substitutions and alterations herein without departing from the scope 
of the invention as defined by the appended claims. 



Claims 

25 

1. A Coriolis mass flow meter, comprising: 

a flow conduit for containing a fluid having a physical characteristic, said fluid adapted to flow in said conduit 
at an unknown rate; 

30 a drive circuit for creating a vibration in said flow conduit, said fluid altering said vibration as a function of said 

physical characteristic and said flow rate; 

a detector circuit for measuring said altered vibration at a working point along a response curve with respect 
to said physical characteristic and producing a signal representing an uncompensated mass flow rate of said 
fluid; and 

35 a computation circuit arranged to calculate a compensated mass flow rate of said fluid based on said uncom- 

pensated rate divided by Cl n , where Q is a natural frequency of said flow conduit at which the conduit is driven, 
and wherein n is an integer or noninteger number except n = -1 , 0 or 1 , chosen in dependence on said working 
point and on the type of sensors used in the detector circuit, said compensated rate thereby being compensated 
for the effects of said physical characteristic. 

40 

2. The meter as recited in Claim 1 wherein said physical characteristic is selected from the group consisting of: 

pressure, 
density and 
4 $ temperature. 

3. The meter as recited in Claim 1 or 2 wherein said working point is selected from the group consisting of: 

time delay, 
so displacement, 
velocity and 
acceleration. 

4. A Coriolis mass flow meter according to claim 1 , wherein the fluid has a pressure and density, said flow conduit 
55 has a predetermined length, radius and wall thickness and yields a response curve with respect to pressure or 

density that intersects zero at a thereby defined working point, thereby rendering said flow conduit substantially 
insensitive to pressure or density at said working point, and 
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the detector circuit for measuring said altered vibration at said working point produces a signal representing 
a partially compensated mass flow rate that, by virtue of said intersection, is rendered substantially independent 
of effects of pressure or density. 

5 5. The meter as recited in Claim 1 or 4 wherein said compensated rate is related to said uncompensated rate by: 

M2 = M1 * B * (A/a n ) 

10 where: 

M2 = said compensated rate; 
M1 = said at least partially uncompensated rate; 
B = temperature and strain compensation values; 
is A = a scaling factor; 

CI = said natural frequency; and 
n = said number. 

6. A Coriolis mass flow meter according to claim 4 wherein said flow conduit yields pressure and density response 
20 curves that both intersect zero at a thereby defined working point, thereby rendering said flow conduit substantially 

insensitive to said pressure and density at said working point; and 

the detector circuit produces a signal representing a mass flow rate that, by virtue of said intersection, is 
rendered substantially independent of effects of said pressure and density. 

25 

7. The meter according to Claim 4 or 6 wherein said working point is along a curve of 1/n n , where ft is a driven 
natural frequency of said flow conduit and n is a number determining said working point and selected from the 
group consisting of: 

30 a positive integer real value, 

a positive noninteger real value, 
a negative integer real value and 
a negative noninteger real value. 

35 8. The meter according to Claim 4,6 or 7 wherein said working point is an integer real value, said detector circuit 
thereby measuring a characteristic of said altered vibration, said characteristic selected from the group consisting 
of: 

time delay, 
^0 displacement, 
velocity and 
acceleration. 

9. The meter according to any of Claims 1 to 8 wherein said drive circuit creates a vibration selected from the group 
45 consisting of: 

a bending mode vibration and 
a radial mode vibration. 

so 10. The meter according to any of Claims 1 to 9 further comprising a temperature sensor coupled to said conduit for 
measuring a temperature thereof and producing a temperature signal related thereto, said computation circuit 
employing said temperature signal to calculate said compensated rate. 

1 1 . The meter-according to any of Claims 1 to 1 0 further comprising a strain gage coupled to said conduit for measuring 
55 a strain thereof and producing a strain signal related thereto, said computation circuit employing said strain signal 

to calculate said compensated rate. 

12. The meter according to any of Claims 1 toll further comprising a pressure-tight case surrounding said flow conduit 
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and defining a space about said conduit for containing a prescribed amount of pressure. 

13. The meter according to any of Claims 1 to 12 further comprising manifolds fixedly attached to ends of said flow 
conduit to allow said fluid to flow into and out of said flow conduit. 

5 

14. The meter according to any of Claims 1 to 1 3 wherein said drive circuit comprises a plurality of magnet/coil pairs 
located along a length of said flow conduit. 

15. A method of determining a mass flow rate of a fluid using a Coriolis mass flow meter, comprising the steps of: 

containing said fluid in a flow conduit, said fluid having a physical characteristic, said fluid adapted to flow in 
said conduit at an unknown rate; 

creating a vibration in said flow conduit with a drive circuit, said fluid altering said vibration as a function of 
said physical characteristic and said flow rate; 

measuring said altered vibration at a working point along a response curve with respect to. said physical 
characteristic with a detector circuit, said detector circuit producing a signal representing an uncompensated 
mass flow rate of said fluid; and 

calculating a compensated mass flow rate ot said fluid with a computation circuit, said compensated rate being 
based on said uncompensated rate divided by Q n , where Q is a natural frequency of said flow conduit at which 
the conduit is driven and wherein n is an integer or noninteger number except n = -1 , 0, or 1, chosen in 
dependence on said working point and on the type of sensors used in the detector circuit, said compensated 
rate thereby being compensated for the effects of said physical characteristic. 

16. The method as recited in Claim 1 5 wherein said step of containing comprises the step of containing a fluid having 
25 a fluid physical characteristic selected from the group consisting of: 

pressure, 
density and 
temperature. 

30 

17. The method as recited in Claim 15 or 16 wherein said step of measuring comprises the step of measuring said 
altered vibration at a working point selected from the group consisting of: 

time delay, 
35 displacement, 
velocity and 
acceleration. 

18. A method of determining a mass flow rate according to claim 15 wherein said fluid has a pressure and density, 
40 said flow conduit has a predetermined length, radius and wall thicknes and yields a response curve that intersects 

zero at a therby defined working point along said response curve with respect to pressure or density, thereby 
rendering said flow conduit substantially insensitive to pressure or density at said working point, and 

said detector circuit produces a signal representing a partially compensated mass flow rate that, by virtue of 
45 said intersection, is rendered substantially independent of effects of pressure or density. 

19. The method as recited in Claim 15 or 18 wherein said step of calculating comprises the step of relating said 
compensated rate to said uncompensated rate by: 

50 

M2 = M1 *B*(A/H n ) 

where: 

55 M2 = said compensated rate; 

M1 = said at least partially uncompensated rate; 
B = temperature and strain compensation values; 
A = a scaling factor; 
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ft = said natural frequency; and 
n = said number. 

20. A method of determining a mass flow rate according to Claim 19 wherein said flow conduit yields pressure and 
5 density response curves that both intersect zero at a thereby defined working point, thereby rendering said flow 

conduit substantially Insensitive to pressure and density at said working point, and 

said detector circuit produces a signal representing a mass flow rate that, by virtue of said intersection, is 
rendered substantially independent of effects of said pressure and density 

10 

21. The method as recited in Claim 18 or 20 wherein said step of creating comprises the step of creating a vibration 
in said flow conduit wherein said working point is along a curve of 1/ft n , where Q. is a driven natural frequency of 
said flow conduit and n is a number determining said working point and selected from the group consisting of: 

is a positive integer real value, 

a positive noninteger real value, 
a negative integer real value and 
a negative noninteger real value. 

20 22. The method according to Claim 1 8,20 or 21 wherein said working point is an integer real value, said step of meas- 
uring comprising the step of measuring a characteristic of said altered vibration with said detector circuit, said 
characteristic selected from the group consisting of: 

time delay, 
25 displacement, 
velocity and 
acceleration. 

23. The method according to any of Claims 15 to 22 wherein said step of creating comprises the step of creating a 
30 vibration selected from the group consisting of: 

a bending mode vibration and 
a radial mode vibration. 

35 24. The method according to any of Claims 15 to 23 further comprising the step of coupling a temperature sensor to 
said conduit for measuring a temperature thereof, said temperature sensor producing a temperature signal related 
thereto, a computation circuit of said meter employing said temperature signal to calculate said compensated rate. 

25. The method according to any of Claims 15 to 24 further comprising the step of coupling a strain gage to said 
40 conduit for measuring a strain thereof, said strain gage producing a strain signal related thereto, said computation 

circuit employing said strain signal to calculate said compensated rate. 

26. The method according to any of Claims 1 5 to 25 further comprising the step of surrounding said flow conduit with 
a pressure-tight case, said case defining a space about said conduit for containing a prescribed amount of pressure. 

45 

27. The method according to any of Claims 1 5 to 26 further comprising the step of fixedly attaching manifolds to ends 
of said flow conduit to allow said fluid to flow into and out of said flow conduit. 

28. The method according to any of Claims 15 to 27 further comprising the step of locating a plurality of magnet/coil 
50 pairs of said drive circuit along a length of said flow conduit. 

29. A method according to Claim 14 wherein the flow conduit contains a fluid having a pressure and a density and is 
composed of a desired material and has a desired length, radius and wall thickness; 

55 simultaneously two natural modes of vibration are created in said flow conduit with the drive circuit, said fluid 

altering said modes of vibration as a function of said pressure, density, flow rate, material, length, radius and 
wall thickness; and 

said detector circuit produces a first signal representing an uncompensated mass flow rate of said fluid; further 
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comprising the steps of: 

changing said pressure and density of said fluid; 

again measuring said altered modes of vibration with the detector circuit at said working point, said detector 
circuit producing a second signal representing an uncompensated mass flow rate of said fluid; 
s computing said pressure and density by simultaneous solution, said simultaneous solution further providing 

a compensation factor relating to said material, length, radius and wall thickness of said flow conduit; and 
using a computation circuit for said meter embodying said compensation factor, to provide a compensated 
mass flow signal. 

10 30. The method as recited in Claim 29 wherein said step of computing comprises the step of simultaneously solving: 

Q = F1 Q + (C1)*P - (C2)*D - (C5)T 

is and 

q = F2 0 + (C3)*P - (C4)*D- (C7)T 

20 where: 

Q and Q - measured frequencies of said first and second signals; 

F1 0 and F2 0 = zero pressure/density frequency values for said modes of vibration, respectively; 
C1 - C10 = predetermined constants; 

25 p = said pressure; and 

D = said density. 

31. The method as recited in Claim 29 or 30 wherein said step of measuring comprises the step of measuring a 
temperature and axial stress of said flow conduit and a viscosity of said fluid. 



30 



3S 



Paten tanspruche 

1. Coriolis-MassendurchfluB-MeBeinrichtung, mit: 



einer DurchfluB-Leitung zur Aufnahme eines Fluids mit einer physikalischen Eigenschaft, wobei das Fluid 
geeignet ist, in der Leitung mit einer unbekannten Rate zu flieBen; 

einer Antriebsschaftung zum Erzeugen einer Vibration in der DurchfluB-Leitung, wobei das Fluid die vibration 
als Funktion der physikalischen Eigenschaft und der DurchfluBrate verandert; 
40 einer Detektor-Schaltung zum Messen der veranderten Vibration bei einem Arbeitspunkt entlang einer die 

physikalische Eigenschaft betreffenden Ansprechkurve und zum Erzeugen eines Signals, das eine unkom- 
pensierte MassendurchfluBrate des Fluids darstellt; und 

einer Berechnungs-Schaltung, die zur Berechnung einer kompensierten MassendurchfluBrate des Fluids ba- 
sierend auf der un kompensierten Rate dividiert durch Q n ausgebildet ist, wobei Q eine Eigenfrequenz der 
45 DurchfluB-Leitung ist, mit der die Leitung angetrieben wird, und wobei n eine ganze oder nicht-ganze Zahl, 

auBer n = -1, 0 oder 1, ist, die abhangig von dem Arbeitspunkt und von dem in der Detektor-Schaltung ver- 
wendeten Sensortyp ausgewahlt wird, wobei dadurch die Wirkungen der physikalischen Eigenschaft in der 
kompensierten Rate kompensiert sind. 

so 2. MeBeinrichtung nach Anspruch 1 , wobei die physikalische Eigenschaft aus der Gruppe ausgewahlt wird, die be- 
steht aus: 

Druck, 
Dichte und 
ss Temperatur. 

3. MeBeinrichtung nach Anspruch 1 oder 2, wobei der Arbeitspunkt aus der Gruppe ausgewahlt wird, die besteht aus: 
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Zeitverzogerung, 
Auslenkung, 
Geschwindigkeit und 
Beschleunigung. 

Coriolis-MassendurchfluG-MeBeinrichtung nach Anspruch 1, wobei das Fluid einen Druck und eine Dichte besitzt, 
die DurchfluG-Leitung eine vorbestimmte Lange, einen Radius und eine Wanddicke aufweist und zu einer Druck 
und Dichte betreff enen Ansprechkurve fflhrt, die den Nullwert an einem dadurch definierten Arbeitspunkt schneidet, 
wobei die DurchfluG-Leitung dadurch an diesem Arbeitspunkt im wesentlichen unempfindlich gegenuber Druck 
und Dichte gemacht wird, und 

die Detektor-Schaltung zum Messen der veranderten Vibration an dem Arbeitspunkt ein Signal erzeugt, das 
eine teilweise kompensierte Massendurchflu Grate darstellt, die aufgrund des Schnittpunkts im wesentlichen 
unabhangig von Wirkungen des Drucks und der Dichte ist. 

MeGeinrichtung nach Anspruch 1 Oder 4, wobei die kompensierte Rate in Bezug zu der unkompensierten Rate 
gesetzt wird durch: 

M2 = M1 * B * (A/Q n ) 

wobei: 

M2 = die kompensierte Rate; 

M1 = die zumindest teilweise kompensierte Rate; 

B = Temperatur und Dehnungskompensationswerte; 

A = ein Skalierungsfaktor; 

CI = die Eigenfrequenz; und 

n = die Zahl. 

Coriolis-MassendurchfluG-MeGeinrichtung nach Anspruch 4, wobei die DurchfluG-Leitung Druck- und Dichte-An- 
sprechkurven aufweist, die beide den Nullpunkt an einem dadurch definierten Arbeitspunkt schneiden, wodurch 
die DurchfluG-Leitung im wesentlichen unempfindlich gegenuber dem Druck und der Dichte an diesem Arbeits- 
punkt wird; und 

die Detektor-Schaltung ein Signal erzeugt, das eine Massendurchflu Grate darstellt, die aufgrund des Schnitt- 
punktes im wesentlichen unabhangig von den Wirkungen des Drucks und der Dichte ist. 

MeGeinrichtung nach Anspruch 4 oder 6, wobei der Arbeitspunkt entlang einer Kurve von 1/Q n liegt, wobei Q. die 
angetriebene Eigenfrequenz der DurchfluG-Leitung ist und n eine Zahl ist, die den Arbeitspunkt bestimmt und aus 
der Gruppe ausgewahlt wird, die umfaGt: 

einen posrtiven ganzzahligen realen Wert, 
einen positiven nicht-ganzzahligen realen Wert, 
einen negativen ganzzahligen realen Wert und 
einen negativen nicht-ganzzahligen realen Wert. 

MeGeinrichtung nach Anspruch 4, 6, 7, wobei der Arbeitspunkt ein ganzzahliger realer Wert ist, wobei die Detektor- 
Schaltung dadurch eine Eigenschaft der veranderten Vibration miGt, wobei die Eigenschaft aus der Gruppe aus- 
gewahlt wird, die besteht aus: 

Zeitverzogerung, 
Auslenkung, 
Geschwindigkeit und 
Beschleunigung. 

MeGeinrichtung nach einem der Anspruche 1 bis 8, wobei die Antriebsschaltung eine Vibration erzeugt, die aus 
der Gruppe ausgewahlt wird, die besteht aus: 
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Biegemodusvibration und 
Radialmodusvibration . 

10. MeGeinrichtung nach einem der Anspruche 1 bis 9, des weiteren mit einem Temperatursensor, der mit der Leitung 
5 zum Messen deren Temperatur gekoppelt ist und zum Erzeugen eines diese betreffenden Temperatursignals, 

wobei die Berechnungs-Schaltung das Temperatursignal zur Berechnung der kompensierten Rate benutzt. 

1 1 . MeGeinrichtung nach einem der Anspruche 1 bis 1 0, des weiteren mit einem Dehnungsmesser, der mit der Leitung 
zum Messen deren Dehnung gekoppelt ist und zum Erzeugen eines diese betreffenden Dehnungssignals, wobei 

io die Berechnungsschattung das Dehnungssignal zur Berechnung der kompensierten Rate einsetzt. 

12. MeGeinrichtung nach einem der Anspruche 1 bis 11, ferner mit einem druckdichten Gehause, das die DurchfiuG- 
Leitung umgibt und das urn die Leitung einen Zwischenraum zur Aufnahme eines vorbestimmten Drucks bildet. 

*5 13. MeGeinrichtung nach einem der Anspruche 1 bis 12, ferner mit Leitungsrohren, die test mit den Enden der Durch- 
fluB-Leitung verbunden sind, um zu ermoglichen, daG das Fluid in die DurchfluG-Leitung hinein und aus dieser 
herausflieGt 

14. MeGeinrichtung nach einem der Anspruche 1 bis 1 3, wobei die Antriebsschaltung eine Vielzahl von Magnet/Spulen- 
20 Paaren umfaGt, die entlang der DurchfluG-Leitung angeordnet sind. 

15. Verfahren zur Bestimmung einer MassendurchfluGrate eines Fluids unter Verwendung einer Coriolis-Massen- 
durchfluG-MeGeinrichtung, mit den Schritten: 

25 Aufnehmen der Flussigkeit in einer DurchfluG-Leitung, wobei das Fluid eine physikalische Eigenschaft besitzt, 

wobei das Fluid geeignet ist, die Leitung mit einer unbekannten Rate zu durchflieGen; 
Erzeugen einer Vibration in der DurchfluG-Leitung mit einer Antriebsschaltung, wobei das Fluid die Vibration 
als Funktion der physikalischen Eigenschaft und der DurchfluGrate verandert; 

Messen der veranderten Vibration an einem Arbeitspunkt entlang einer die physikalische Eigenschaft betref- 
30 f enden Ansprechkurve mit einer Detektor-Schaitung, wobei die Detektor-Schaltung ein Signal erzeugt, das 

eine unkompensierte MassendurchfluGrate des Fluids darstellt; und 

Berechnen einer kompensierten MassendurchfluGrate des Fluids mit einer Berechnungs-Schaltung, wobei 
die kompensierte Rate auf der unkompensierten Rate dividiert durch £l n basiert, wobei Q. die Eigenfrequenz 
der DurchfluG-Leitung darstellt, mit der die Leitung angetrieben wird, und wobei n eine ganze Zahl Oder eine 
3S nicht-ganze Zahl, auGer n = -1, 0, 1, ist und in Abhangigkeit von dem Arbeitspunkt und dem in der Detektor- 

Schaltung eingesetzten Sensortyp ausgewahlt wird, wobei die Wirkung der physikalischen Eigenschaft in der 
kompensierten Rate dadurch kompensiert wird. 

16. Verfahren nach Anspruch 15, wobei der Schritt des Aufnehmen den Schritt des Aufnehmens eines Fluids umfaGt, 
40 das eine physikalische Fluideigenschaft besitzt, die aus der Gruppe ausgewahlt wird, die besteht aus: 

Druck, 
Dichte und 
Temperatur. 

45 

17. Verfahren nach Anspruch 15 oder 16, wobei der Schritt des Messens den Schritt des Messens der veranderten 
Vibration an einem Arbeitspunkt umfaGt, der aus der Gruppe ausgewahlt wird, die besteht aus: 

Zeitverzogerung, 
so Auslenkung, 

Geschwindigkeit und 
Beschleunigung. 

18. Verfahren zum Bestimmen einer MassendurchfluGrate nach Anspruch 15, wobei das Fluid einen Druck und eine 
£5 Dichte aufweist, die DurchfluG-Leitung eine vorbestimmte Lange, einen Radius und eine Wanddicke besitzt und 

eine Ansprechkurve hat, die den Nullpunkt an einem dadurch definierten Arbeitspunkt entlang der Ansprechkurve 
bezuglich des Drucks oder der Dichte schneidet, wodurch die DurchfluG-Leitung im wesentlichen unempfindlich 
gegenuber Druck und Dichte an diesem Arbeitspunkt gemacht wird, und die Detektor-Schaltung ein Signal erzeugt, 
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das eine teilweise kompensierte MassendurchfluGrate darstellt, die aufgrund des Schnittpunkts im wesentlichen 
unabhangig von den Wirkungen des Drucks und der Dichte gemacht wird. 

19. Verfahren nach Anspruch 15 Oder 18, wobei der Schritt des Berechnens den Schritt des in Beziehung Setzens 
der kompensierten Rate zu der unkompensierten Rate umfaGt, durch: 



M2 = M1 *B*(A/Q n ) 



10 wobei 



M2 = die kompensierte Rate; 
M1 = die zumindest teilweise un kompensierte Rate; 
B = Temperatur und Dehnungskompensationswerte; 
15 A = ein Skalisierungsfaktor; 

ft = die Eigenfrequenz; und 
n = die Zahl. 

20. Verfahren zur Bestimmung einer MassendurchfluGrate nach Anspruch 19, wobei die DurchfluG-Leitung zu Druck- 
20 und Dichte-Ansprechkurven f uhrt, die beide den Nullpunkt an einem dadurch definierten Arbeitspunkt schneiden, 

wodurch die DurchfluB-Leitung an dem Arbeitspunkt im wesentlichen unempfindlich gegenuber Druck und Dichte 
gemacht wird, und 

die Detektor-Schaltung ein Signal erzeugt, das die MassendurchfluGrate darsteflt, die aufgrund des Schnitt- 
25 punkts im wesentlichen unabhangig von den Wirkungen von Druck und Dichte gemacht wird. 

21. Verfahren nach Anspruch 18 Oder 20, wobei der Schritt des Erzeugens den Schritt des Erzeugens einer Vibration 
in der DurchfluG-Leitung umfaGt, wobei der Arbeitspunkt entlang einer Kurve von 1/Q n liegt, wobei ft eine ange- 
triebene Eigenfrequenz der DurchfluG-Leitung ist und n eine Zahl ist, die den Arbeitspunkt bestimmt und aus der 

30 Gruppe gewahlt wird, die umfaGt: 

einen positiven ganzzahligen realen Wert, 
einen positiven nicht-ganzzahligen realen Wert, 
einen negativen ganzzahligen realen Wert, und 
35 einen negativen nicht-ganzzahligen realen Wert. 

22. Verfahren nach Anspruch 1 8, 20 Oder 21 , wobei der Arbeitspunkt ein ganzzahliger realer Wert ist, wobei der Schritt 
des Messens den Schritt des Messens einer Eigenschaft der veranderten Vibration mit der Detektions-Schaltung 
umfaGt, wobei die Eigenschaft aus der Gruppe ausgewahlt wird, die besteht aus: 

40 

Zeitverzogerung, 
Auslenkung, 
Geschwindigkeit und 
Beschleunigung. 

45 

23. Verfahren nach einem der Anspruche 1 5 bis 22, wobei der Schritt des Erzeugens den Schritt des Erzeugens einer 
Vibration umfaGt, die aus der Gruppe ausgewahlt wird: 

Biegemodusvibration, und 
so Radialmodusvibration. 

24. Verfahren nach einem der Anspruche 15 bis 23, ferner mit dem Schritt des Koppelns eines Temperatursensors 
mit der Leitung zum Messen deren Temperatur, wobei der Temperatursensor ein diese betreffendes Temperatur- 
signal erzeugt, wobei eine Berechnungs-Schaltung der MeGeinrichtung das Temperatursignal zur Berechnung der 

55 kompensierten Rate benutzt. 

25. Verfahren nach einem der Anspruche 1 5 bis 24, ferner mit dem Schritt des Koppelns eines Dehnungsmessers mit 
der Leitung zum Messen deren Dehnung, wobei der Dehnungsmesser ein diese betreffendes Dehnungssignal 
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erzeugt, wobei die Berechnungs-Schaltungdas Dehnungssignal zur Berechnungder kompensierten Rate benutzt. 

26. Verfahren nach einem der Anspruche 15 bis 25, femer mit dem Schritt des UmschlieBens der DurchfluB-Leitung 
mit einem druckdichten Gehause, wobei das Gehause urn die Leitung einen Raum zur Aufnahme eines vorgege- 

5 benen Drucks ausbildet. 

27. Verfahren nach einem der Anspruche 15 bis 26, ferner mit dem Schritt der festen Anbringung von Leitungsrohren 
an Enden der DurchfluB-Leitung, um dem Fluid zu ermoglichen, in die DurchfluB-Leitung hinein und heraus zu 
flieBen. 

10 

28. Verfahren nach einem der Anspruche 1 5 bis 27, ferner mit dem Schritt der Anordnung einer vielzahl von Magnet/ 
Spulen-Paaren der Antriebs-Schaltung entlang der DurchfluB-Leitung. 

29. Verfahren nach Anspruch 14, wobei die DurchfluB-Leitung ein Fluid enthalt, das einen Druck und eine Dichte 
f5 aufweist und sich aus einem gewunschten Material zusammensetzt und eine gewunschte Lange, einen Radius 

und eine Wanddicke besitzt; 

wobei gleichzeitig zwei Vibrations-Eigenmodi in der DurchfluB-Leitung mit der Antriebs-Schaltung erzeugt 
werden, wobei das Fluid die Vibration smodi als Funktion des Drucks, der Dichte, der DurchfiuBrate, des Ma- 
20 terials, der Lange, des Radius und der Wanddicke verandert; und 

die Detektor-Schaltung ein erstes Signal erzeugt, das eine unkompensierte MassendurchfluBrate des Fluids 

darstellt; ferner mit den Schritten: 

Andern des Drucks und der Dichte des Fluids; 

erneutes Messen der veranderten Vibrationsmodi mit der Detektor-Schaltung an dem Arbeitspunkt, wobei die 
25 Detektor-Schaltung ein zweites Signal erzeugt, das eine unkompensierte MassendurchfluBrate des Fluids 

darstellt; 

Berechnen des Drucks und der Dichte durch eine gleichzeitige Auflbsung, wobei die gleichzeitige Auflosung 
ferner einen Kompensationsfaktor iiefert, der das Material, die Lange, den Radius und die Wanddicke der 
DurchfluB-Leitung betrifft; und 

so Benutzen einer Berechnungs-Schaltung fur die MeBeinrichtung, die den Kompensationsfaktor bildet, um ein 

kompensiertes MassendurchfluBsignal zu liefern. 

30. Verfahren nach Anspruch 29, wobei der Schritt der Berechnung den Schritt der gleichzeitigen Auflosen umfaBt: 

Q = F1 0 + (C1)*P - (C2)*D - (C5)T 



40 



so 



und 



und wobei: 



a = F2 0 + (C3)*P - (C4)*D- (C7)T 



45 n und Q - gemessene Frequenzen des ersten und des zweiten Signals; 

F1 0 und F2 0 = Nulldruck- bzw. Dichtefrequenzwerte der Vibrationsmodi; 

C1 - C10 = vorbestimmte Konstanten; 

P = der Druck; und 

D = die Dichte. 



31. Verfahren nach Anspruch 29 Oder 30, wobei der Schritt des Messens den Schritt des Messens einer Temperatur 
und einer axialen Spannung der DurchfluB-Leitung und einer Viskositat des Fluids umfaBt. 



55 Revendfcations 

1 . Debitmetre massique a effet de Coriolis, comportant: 
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un conduit d'ecoulement destine a contenir un fluide ayant une caracteristique physique, ledit fluide etant 
prevu pour s'ecouler dans ledit conduit a un debit inconnu; 

un circuit d'excitation destine a produire une vibration dans ledit conduit d'ecoulement, ledit fluide modifiant 
ladite vibration en fonction de ladite caracteristique physique et dudit debit d'ecoulement; 
un circuit de detecteur destine a mesurer ladite vibration modifiee en un point de travail le long d'une courbe 
de reponse par rapport a ladite caracteristique physique et produire un signal representant un debit massique 
non compensg dudit fluide; et 

un circuit de calcul prevu pour calculer un debit massique compense dudit fluide sur la base dudit debit non 
compensd divise par QP, ou Q est une frequence propre dudit conduit d'ecoulement a laquelle le conduit est 
excite, et n est un nombre entier ou un nombre non entier, excepte n = -1 , 0 ou 1 , choisi en fonction dudit point 
de travail et du type de capteurs utilises dans le circuit de detecteur, ledit debit compense etant ainsi compense 
des effets de ladite caracteristique physique. 

Debitmetre selon la revendication 1, dans lequel ladite caracteristique physique est choisie dans le groupe com- 
pose de : 

la pression, 
la densite, et 
la temperature. 

Debitmetre selon la revendication 1 ou 2, dans lequel ledit point de travail est choisi dans le groupe compose de : 

le retard temporel, 
le deplacement, 
la Vitesse, et 
I'accele ration. 

Debitmetre massique a eff et de Coriolis selon la revendication 1 , dans lequel le fluide a une pression et une densite, 
ledit conduit d'ecoulement a une longueur, un rayon et une epaisseurde paroi predetermines et donne une courbe 
de reponse par rapport a la pression ou la densite qui passe par z6ro en un point de travail ainsi defini, rendant 
de ce fait ledit conduit d'ecoulement sensiblement insensible a la pression ou la densite au niveau dudit point de 
travail, et 

le circuit de detecteur destine a mesurer ladite vibration modified au niveau dudit point de travail produit un 
signal representant un debit massique partieilement compense qui, en vertu de ladite intersection, est rendue 
sensiblement independant des effets de la pression ou de la densite. 

Debitmetre selon la revendication 1 ou 4, dans lequel ledit debit compense est lie audit debit non compense par: 

M2 = M1 * B * (A/Q n ) 

oil : 

M2 = ledit debit compense; 

M1 = ledit d6bit au moins partieilement non compens6; 

B = valeurs de compensation de temperature et de contrainte; 

A = un facteur d'echelle; 

& = ladite frequence propre; et 

n = ledit nombre. 

Debitmetre massique a effet de Coriolis selon la revendication 4, dans lequel ledit conduit d'ecoulement donne 
des courbes de teponse de pression et de densite qui passent toutes les deux par zero en un point de travail ainsi 
defini, rendant de ce fait ledit conduit d'ecoulement en grande partie insensible a la ladite pression et a ladite 
densite audit point de travail; et 

le circuit de detecteur produit un signal representant un debit massique qui, en vertu de ladite intersection, 
est rendu en grande partie independant des effets de ladite pression et de ladite densite. 
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7. Debitmetre selon la revendication 4 ou 6, dans lequel ledit point de travail est situ6 le long d'une courbe de 1/Q n , 
ou Q est une frequence propre excitee dudit conduit d'ecoulement et n est un nombre determinant ledit point de 
travail et choisi dans le groupe compose de : 

5 une valeur r6elle entiere positive, 

une valeur reelle non entiere positive, 
une valeur reelle entiere negative, et 
une valeur reelle non entiere negative. 

10 8. Debitmetre selon la revendication 4, 6 ou 7, dans lequel ledit point de travail est une valeur r6elle entiere, ledit 
circuit de detecteur mesurant ainsi une caracteristique de ladite vibration modiftee, ladite caracteristique etant 
choisie parmi le groupe compose de : 

le retard temporel, 
15 le deplacement, 

la Vitesse, et 
{'acceleration. 

9. Debitmetre selon Tune quelconque des revendications 1 a 8, dans lequel ledit circuit d'excitation produit une vi- 
20 brat ion choisie dans le groupe compose de : 

une vibration de mode de flexion, et 
une vibration de mode radial. 

25 10. Debitmetre selon Tune quelconque des revendications 1 a 9, comportant en outre un capteur de temperature relie 
audit conduit afin de mesurer une temperature de celui-ci et produire un signal de temperature qui s'y rapporte, 
ledit circuit de calcul utilisant ledit signal de temperature afin de calculer ledit debit compense. 

11. Debitmetre selon Tune quelconque des revendications 1 a 10, comportant en outre une jauge de contrainte reliee 
30 audit conduit afin de mesurer une contrainte de celui-ci et produisant un signal de contrainte qui s*y rapporte, ledit 

circuit de calcul utilisant ledit signal de contrainte afin de calculer ledit debit compense. 

12. Debitmetre selon Tune quelconque des revendications 1 a 11 , comportant en outre un cottier resistant a la pression 
entourant ledit conduit d'ecoulement et definissant un espace autour dudit conduit afin de contenir une quantite 

35 prescrite de pression. 

13. Debitmetre selon I'une quelconque des revendications 1 a 12, comportant en outre des collecteurs fix6s a des 
extremites dudit conduit d'6coulement afin de permettre audit fluide de s'ecouler dans et hors dudit conduit d'ecou- 
lement. 

40 

14. Debitmetre selon Tune quelconque des revendications 1 a 13, dans lequel ledit circuit d'excitation comporte plu- 
sieurs paires aimant/bobine situ6es sur une longueur dudit conduit d'ecoulement. 

15. Precede de determination d'un debit massique d'un fluide a Paide d'un debitmetre massique a effet de Coriolis, 
45 comportant les etapes consistant a : 

contenir ledit fluide dans un conduit d'ecoulement, ledit fluide ayant une caracteristique physique, ledit fluide 
etant prevu pour s'ecouler dans ledit conduit a un debit inconnu; 

produire une vibration dans ledit conduit d'ecoulement avec un circuit d'excitation, ledit fluide modifiant ladite 
50 vibration en fonction de ladite caracteristique physique et dudit debit d'ecoulement; 

mesurer ladite vibration modifiee en un point de travail le long d'une courbe de reponse par rapport a ladite 
caracteristique physique avec un circuit de detecteur, ledit circuit de detecteur produisant un signal repr6sen- 
tant un debit massique non compense dudit fluide; et 

calculer un debit massique compense dudit fluide avec un circuit de calcul, ledit debit compense etant bas6 
55 sur ledit debit non compense divise par fi n , ou Q est une frequence propre dudit conduit d'ecoulement a 

laquelle le conduit est excite, et n est un nombre entier ou non entier, excepte n = -1, 0, 1, choisi en fonction 
dudit point de travail et du type de capteurs utilises dans le circuit de detecteur, ledit debit compense etant 
ainsi compense des effets de ladite caracteristique physique. 
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16. Procede selon la revendication 1 5, dans lequel ladite etape destinee a contenir comporte I'etape destinee a contenir 
un fluide ayant une caracteristique physique de fluide choisie dans le groups compose de : 

la pression, 
5 la densite, et 

la temperature. 

17. Procede selon la revendication 15 ou 16, dans lequel ladite etape de mesure comporte I'etape de mesure de tadite 
vibration modifiee en un point de travail choisi dans le groupe compose de : 

10 

le retard temporel, 
le deplacement, 
la vitesse, et 
('acceleration. 

15 

18. Procede de determination d'un debit massique selon la revendication 15, dans lequel ledit fluide a une pression 
et une densite, ledit conduit d'ecoulement a une longueur, un rayon et une epaisseur de paroi predetermines et 
donne une courbe de reponse qui passe par zero en un point de travail ainsi defini le long de ladite courbe de 
reponse par rapport a la pression ou la densite, rendant de ce fait ledit conduit d'ecoulement en grande partie 

20 insensible a la pression ou a ia densite au niveau dudit point de travail, et 

ledit circuit de detecteur produit un signal representant un debit massique partiellement compense qui, en 
vertu de ladite intersection, est rendu en grande partie independant des effets de la pression ou de la densite. 

25 19. Procede selon la revendication 15 ou 18, dans lequel ladite etape de calcul comporte I'etape consistant a Her ledit 
debit compense audit debit non compense par : 



M2 = M1 *B*(A/O n ) 

30 

ou : 

M2 = ledit debit compense; 

M1 = ledit deM au moins partiellement non compense; 
35 B = valeurs de compensation de temperature et de contrainte; 

A = un facte ur d'echelle; 
Q. - ladite frequence propre; et 
n = ledit nombre. 

40 20. Procede de determination d'un debit massique selon la revendication 19, dans lequel ledit conduit d'ecoulement 
donne des courbes de reponse de pression et de densite qui passent toutes les deux par zero en un point de 
travail ainsi defini, rendant de ce fait ledit conduit d'ecoulement en grande partie insensible a la pression et a la 
densite au niveau dudit point de travail, et 

45 ledit circuit de detecteur produit un signal representant un debit massique qui, en vertu de ladite intersection, 

est rendu en grande partie independant des effets de ladite pression et de ladite densite. 

21. Procede selon la revendication 18 ou 20, dans lequel ladite etape de production comporte I'etape de production 
d'une vibration dans ledit conduit d'ecoulement, ledit point de travail etant situe le long d'une courbe de 1/Q n , ou 

so a est une frequence propre excitee dudit conduit d'ecoulement et n est un nombre determinant ledit point de travail 

et choisi dans le groupe compose de : 

une valeur reelle entiere positive, 
une valeur reelle non entiere positive, 
55 une valeur reelle entiere negative, et 

une valeur reelle non entiere negative. 

22. Procede selon la revendication 18, 20 ou 21 , dans lequel ledit point de travail est une valeur reelle entiere, ladite 
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etape de mesure comportant I'etape de mesure d'une caracteristique de ladite vibration modifiee avec led it circuit 
de detecteur, ladite caracteristique etant choisie dans le groupe compost de : 

le retard tempore!, 
5 le deplacement, 

la Vitesse, et 
I'acceleration. 

23. Precede selon Tune quelconque des revendications 15 a 22, dans lequel ladite etape de production comporte 
10 I'etape de production d'une vibration choisie dans le groupe compose de : 

une vibration de mode de flexion, et 
une vibration de mode radial. 

is 24. Precede selon Tune quelconque des revendications 15 a 23, comportant en outre I'etape consistant a relier un 
capteur de temperature audit conduit afin de mesurer une temperature de celui-ci, ledit capteur de temperature 
produisant un signal de temperature qui s'y rapporte, un circuit de calcul dudit d6bitmetre utilisant ledit signal de 
temperature pour calculer ledit debit compense. 

20 25. Procede selon I'une quelconque des revendications 15 a 24, comportant en outre I'etape consistant a relier une 
jauge de contrainte audit conduit afin de mesurer une contrainte de celui-ci, ladite jauge de contrainte produisant 
un signal de contrainte qui s'y rapporte, ledit circuit de calcul utilisant ledit signal de contrainte pour calculer ledit 
d6bit compense. 

25 26. Procede selon Tune quelconque des revendications 15 a 25, comportant en outre I'etape consistant a entourer 
ledit conduit d'ecoulement avec un boitier resistant a la pression, ledit boTtier definissant un espace autour dudit 
conduit afin de contenir une quantite prescrite de pression. 

27. Precede selon Tune quelconque des revendications 15 a 26, comportant en outre I'etape consistant a fixer des 
30 collecteurs a des extremites dudit conduit d'ecoulement afin de permettre audit fluide de s'6couler dans et hors 

dudit conduit d'ecoulement. 

28. Procede selon Tune quelconque des revendications 15 a 27, comportant en outre I'etape consistant a placer plu- 
sieurs paires aimant/bobine dudit circuit d'excitation sur une longueur dudit conduit d'ecoulement. 

35 

29. Precede selon la revendication 14, dans lequel le conduit d'ecoulement contient un fluide ayant une pression et 
une densite et se compose d'une matiere souhaitee et a une longueur, un rayon et une epaisseur de paroi 
souhaites ; 

40 deux modes de vibration propres sont produits simultanement dans ledit conduit d'ecoulement avec le circuit 

d'excitation, ledit fluide modifiant lesdits modes de vibration en fonction de ladite pression, de ladite densite, 
dudit debit, de ladite matiere, de ladite longueur, dudit rayon et de ladite epaisseur de paroi; et 
ledit circuit de detecteur produit un premier signal repr6sentant un debit massique non compense dudit fluide; 

45 comportant en outre les etapes consistant a : 

modifier ladite pression et iadite densite dudit fluide; 

mesurer a nouveau lesdits modes modifies de vibration avec le circuit de detecteur audit point de travail, ledit 
circuit de detecteur produisant un deuxieme signal repr6sentant un debit massique non compense dudit fluide; 
50 calculer ladite pression et ladite densite par resolution simultanee, ladite resolution simultanee procurant en 

outre un facteurde compensation se rapportant a ladite matiere, ladite longueur, ledit rayon et ladite epaisseur 
de paroi dudit conduit d'ecoulement; et 

utiliser un circuit de calcul pour ledit d6bitmetre mettant en oeuvre ledit facteur de compensation, afin de 
procurer un signal de debit massique compense. 

55 

30. Precede selon la revendication 29, dans lequel ladite etape de calcul comporte I'etape de resolution simultanee ; 
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55 



Q = F1 0 + (C1)*P - (C2)*D - (C5)T 

et 

Q = F2 0 + (C3)*P - (C4)*D- (C7)T 

ou 



Q et O = frequences mesur£es desdits premier et deuxieme signaux; 

F1 0 et F2 0 = valeurs de frequence de pression/densitd nulles pour lesdits modes de vibration, respectivement; 
C1 - C10 = constantes pr6d6terminees; 
P = ladite press ion; et 
is D = ladite density. 

31. Procede selon la revendication 29 ou 30, dans lequel ladite etape de mesure comporte l'6tape de mesure d'une 
temperature et d'une contrainte axiale dudit conduit d'ecoulement et d'une viscosit6 dudit fluide. 

20 
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